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Abstract 

'-£3 We present and apply a novel scheme for studying photon- and pion-nucleon scatter- 

ing beyond the threshold region. Partial-wave amplitudes for the 7 TV" and it N states 
are obtained by an analytic extrapolation of subthreshold reaction amplitudes com- 
puted in chiral perturbation theory, where the constraints set by electromagnetic- 
gauge invariance, causality and unitarity are used to stabilize the extrapolation. 

£S| Based on the chiral Lagrangian we recover the empirical s- and p-wave amplitudes 

up to energies y/s ~ 1300 MeV in terms of the parameters relevant at order Q 3 . 
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1 Introduction 

The study of photon- and pion-nucleon interactions has a long history in 
hadron physics. In recent years such reactions have been successfully used 
as a quantitative challenge of chiral perturbation theory (%PT), which is a 
systematic tool to learn about low-energy QCD dynamics [1..2.3J. The appli- 
cation of %PT is limited to the near threshold region. The pion-nucleon phase 
shifts have been analyzed in great depth at subleading orders in the chiral 
expansion [4,5,6]. Pion photoproduction was studied in [7f8f9]rT0] . Compton 
scattering was considered in [HIT]. 

Though there are many model computations (see e.g. |f 2|ll3|ll4l)15pi6fl7l)18j ) 
that address such reactions at energies significantly larger than threshold it 
is an open challenge to further develop systematic effective field theories that 
have a larger applicability domain than xPT and are predictive neverthe- 
less. The inclusion of the isobar as an explicit degree of freedom in the chiral 
Lagrangian is investigated in [19.20.21.22,23,24,25,26,27]. This leads to an ex- 
tension of the applicability domain of the chiral Lagrangian which is based on 
power counting rules. A successful description of scattering data in the iso- 
bar region requires the systematic summation of an infinite number of terms. 
Such a summation may be motivated by generalized counting rules [23] to 
be applied directly to the S-matrix. Alternatively, the required summation 
may be justified by the request that the process is described in accordance 
with the unitarity constraint [22, 28, 29, 30, 3 lj. The counting rules are applied 
to irreducible diagrams only. An infinite number of reducible diagrams be- 
ing summed by the unitarity request [22J. This is analogous to the scheme 
proposed by Weinberg for the nucleon-nucleon scattering problem [3"2~f3"3"] . 

The purpose of this work is to develop a unified description of photon and 
pion scattering off the nucleon based on the chiral Lagrangian. We aim at 
a description from threshold up to and beyond the isobar region in terms 
of partial-wave amplitudes that are consistent with the constraints set by 
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causality and unitarity. Our analysis is based on the chiral Lagrangian with 
pion and nucleon fields truncated at order Q 3 . We do not consider an explicit 
isobar field in the chiral Lagrangian. The physics of the isobar resonance enters 
our scheme by an infinite summation of higher order counter terms in the 
chiral Lagrangian. The particular summation is performed in accordance with 
unitarity and causality. 

Our work is based on a scheme proposed in [34]. We develop a suitable ex- 
tension to be applied to pion-nucleon scattering, pion photoproduction and 
Compton scattering. The scheme is based on an analytic extrapolation of 
subthreshold scattering amplitudes that is controlled by constraints set by 
electromagnetic-gauge invariance, causality and unitarity. Unitarized scatter- 
ing amplitudes are obtained which have left-hand cut structures in accordance 
with causality. The latter are solutions of non-linear integral equations that are 
solved by N / D techniques. The integral equations are imposed on partial- wave 
amplitudes that are free of kinematical zeros and singularities. Such ampli- 
tudes are constructed in the Appendix. An essential ingredient of the scheme 
is the analytic continuation of the generalized potentials that determine the 
partial-wave amplitudes via the non-linear integral equation. We discuss the 
analytic structure of the generalized potentials in detail and construct suitable 
conformal mappings in terms of which the analytic continuation is performed 
systematically. Contributions from far distant left-hand cut structures are rep- 
resented by power series in the conformal variables. 

The relevant counter terms of the Lagrangian are adjusted to the empirical 
data available for photon and pion scattering off the nucleon. We focus on the 
s- and p-wave partial-wave amplitudes and do not consider inelastic channels 
with two or more pions. We recover the empirical s- and p-wave pion-nucleon 
phase shifts up to about 1300 MeV quantitatively. The pion photoproduction 
process is analyzed in terms of its multipole decomposition. Given the sig- 
nificant ambiguities in those multipoles we offer a more direct comparison of 
our results with differential cross sections and polarization data. A quantita- 
tive reproduction of the data set up to energies of about ~ 1300 MeV is 
achieved. 
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2 Analytic extrapolation of subthreshold scattering amplitudes 



In this section we perform a systematic analytic continuation of the reactions 
ttN — > 7 AT, ttN and 7 A" — > 7 AT, ttN based on the chiral Lagrangian [5f2] . After 
a review of the interaction terms relevant at the accuracy level Q 3 we construct 
the tree-level amplitudes for the three reactions. This specifies our notations 
and conventions. Explicit results for the one-loop diagrams that complement 
the tree-level expressions to chiral order Q 3 are recalled in Appendix A. In 
section 2.1 a non-linear integral equation is presented that combines the con- 
straints of unitarity and causality in terms of a generalized potential. Solutions 
thereof will be constructed in application of the N/D technique. Possible am- 
biguities of N/D technique and their relations to bare resonances or ghosts are 
discussed in section 2.2. In section 2.3 we derive a faithful approximation of the 
generalized potential by means of suitably constructed conformal mappings. 
Electromagnetic gauge invariance is kept rigourously. 

The starting point of our studies is the chiral Lagrangian involving pion, nu- 
cleon and photon fields [5f2] . For the readers' convenience we collect all terms 
that contribute to order Q 3 for pion-nucleon scattering, pion photoproduction 
and Compton scattering. From [5]fTT] we extract the terms 



Ant = -^p N>f (r • (7? x (0 M #))) N + |y AT 75 7^ (f • (d^)) N 

- e { (7? x («Vf)) 3 + i±p N - iV 75 7^ (f x tt) 3 iv} A" 

Am N 2 32vr 2 / 3 ^ 

- ^ ml N (7? • 7?) AT - — (iV (0„ 7?) ■ (3„t?) (d^d^N) + h.c. 

+ ^ iV (3„ 7?) • (c^t?) N-^Na^(f- ((^tt) x (9,7?))) iV 

- i iV (r ■ ((^7?) x (9^7?))) (d"N) + h.c. 

+ iV (r ■ ((^7?) x (9^7?))) {d"d»d x N) + h.c. 

2 j 

- 2 z A" (r • (7? x (^7?))) (d^N) + h.c. 

- e^ a/? A 7 (d 8 (0 a vr 3 ) + rf 9 (f • (9 a 7?))) (fy iV) F^ + h.c. 
+ 1 d U ~ dl5 Na^ Ud„nf) ■ (9^ A 7?)) (9 A iV) + h.c. 
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Table 1 

Isospin coefficients as introduced in Q. 
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with the field strength tensor F^ v = d^A v — d u A^. We use the convention 
e 0123 = +1. Most of the parameters introduced in ([!]) are reasonably known 
from xPT studies of the pion-nucleon scattering and pion-photoproduction 
processes close to threshold. We will connect to the different parameter sets 
suggested in the literature in the result section when discussing our preferred 
parameter set. 

In a strict chiral expansion the order Q 3 results are composed from a tree- 
level part and a one-loop part. Both parts are invariant under electromagnetic 
gauge transformations separately. 

The tree-level pion-nucleon scattering amplitude receives contributions from 
the Weinberg- Tomozawa term, the s- and u-channel nucleon exchange pro- 
cesses and the Q 2 and Q 3 counter terms characterized by the parameters 
ci, C4 and di + d 2 , d$, d$, du — d± 5 , dig. We identify the relevant terms 
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Table 2 

Isospin coefficients as introduced in Q. 
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with the nucleon propagators Sn(p), defined as follows 



= I 2 1 ■ • 3 

We use p^ and e/ M (q^) for the initial (final) 4-momenta of the baryon 
and meson. The pion-nucleon-scattering amplitudes are projected onto good 
isospin I where the various isospin coefficients are collected in Tab. [I] 
The contribution from the one-loop diagrams is recalled in Appendix A. 

At tree-level the pion photoproduction amplitude is determined by the Kroll- 
Rudermann term, the nucleon s- and u-channel exchange processes and the 
t-channel pion exchange. At chiral order Q 3 the counter terms d 8 ,d Q , di 8 , 2 d 2 \ — 
d 2 2 contribute. It holds 
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The photoproduction amplitude depends on the choice of the nucleon (p, n) of 
the initial state: the amplitude t!£^n * s obtained from T^^^- by replacing 
the corresponding isospin coefficients of Tab. [2] 



The proton Compton scattering is given by nucleon s- and u-channel exchange 
contribution and pion t-channel exchange. 
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with the proton's anomalous magnetic moment k p = Ka + Kv . 
2. 1 Constraints from unitarity and causality 



Our approach is based on partial-wave dispersion relations, for which the uni- 
tarity and causality constraints can be combined in an efficient manner. Using 
a partial-wave decomposition simplifies calculations because of angular mo- 
mentum and parity conservation. For a suitably chosen partial-wave ampli- 
tude with angular momentum J, parity P and channel quantum numbers a, b 
(in our case ttN and 7iV channels) we separate the right-hand cuts from the 
left-hand cuts 



ifth, TT W - fiM W - y/S - 16 

where the generalized potential, U^ P \y/s), contains left-hand cuts only, by 
definition. It is emphasized that the separation (|6]) is gauge invariant. This fol- 
lows since both contributions in ^ are strictly on-shell and characterized by 
distinct analytic properties]^] The amplitude is considered as a function of s/s 



1 Gauge invariant approximations for the Bethe-Salpeter equation are much more 
difficult to construct |35j. This is because the scattering kernel is required for off- 
shell kinematics and therefore is gauge dependent in this case. 
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due to the MacDowell relations [36J. A subtraction is made at y/s = pu — ttin 
for the reasons to be discussed below. The relation ^ illustrates that the am- 
plitude possesses a unitarity cut along the positive real axis starting from the 
lowest s-channel threshold. In our case the jN intermediate states induces a 
branch point at y/s = p t hr — mN, which defines the lowest s-channel unitarity 
threshold. The structure of the left-hand cuts in U^ P \y/s) can be obtained by 
assuming the Mandelstam representation [37f38] or examining the structure 
of Feynman diagrams in perturbation theory [39] . The form of the analyticity 
domain of partial-wave amplitudes as implied by basic principles of quantum 
field theory is discussed in (ID] . Without specifying a particular structure of the 
left-hand singularities one can regard (6) definition of U ( J P \y/s) |2TK2] . 
In what follows we will assume for definiteness the Mandelstam representation 
to hold. 

The condition that the scattering amplitude must be unitary allows one to 
calculate the discontinuity along the right hand cut 



ATjf'tVS) = 1 (Tjf'tVS + U) - IT'IVS - 

= E ^(v^ + ^plf \Vs)T% P \V~s-ie), (7) 



c,d 



where, p^ 1 (y/s), is the phase-space matrix. The sum in (7) runs over all 
possible intermediate states. Combining ^ and ^ we arrive at the following 
non-linear integral equation for the coupled-channel scattering amplitude 



rp / rz \ tt 1 /tn , V- [ dw Vs- p M T ac (w) p cd (w) T db (w) 

7^ ^Mthra TT W - p M W - y/S -16 



where we suppressed reference to angular momentum and parity. Given any 
gauge- invariant approximation to the generalized potential U a b(y/s), the non- 
linear integral equation ^ generates a gauge-invariant coupled-channel scat- 
tering amplitude. 

It is convenient to use partial-wave amplitudes free of kinematic constraints, 
such as zeros and singularities at threshold and pseudothreshold. In Ap- 
pendix [A] partial-wave amplitudes are established, that are kinematically un- 
constrained. The partial-wave amplitudes are scaled by powers of yfs such 
that the phase- space factors p a b(y/s) approach finite values for asymptoti- 
cally large energies. This implies that the partial-wave scattering amplitudes 
are bounded asymptotically in all partial waves. The generalized potentials 
are bounded modulo some logarithmic structures. All together we derive 
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with the total angular momentum J and parity P 

for J — | even. In (9D we consider p cm and £7, the relative momentum and the 
nucleon energy in the center-of-mass system, as a function of Mandelstam's 
variable s. 



Given an approximate generalized potential we will solve for the partial-wave 
scattering amplitudes ^ with the phase-space distributions In the present 
rather exploratory study we take into account only the 7rN intermediate states 
in ([8]). The higher mass intermediate states are then effectively included in 
the potential uj£ P \y/s). We also neglect the 'yN intermediate state, since the 
effects thereof are suppressed by the square of the electromagnetic charge. 

The non-linear integral equation ^ can be solved by means of N / D techniques 
|13] . This problem was studied extensively in the literature j4*4l45l!46l)47!l42] . 
The amplitude is represented as 



T ab (V~S ) = £ D ac\V~S ) AUV^ ) , (10) 

c 

where D ab (y/s) has no singularities but the right-hand s-channel unitarity 
cuts. The coupled-channel unitarity condition is a consequence of the ansatz 



n ( r\ x ^ f°° dw - fi M N ac (w) p cb (w) 
^ , ab{\ / s) = O ab -y j= — . (11) 

V-Vthr 7T W - II M W- ' 



In contrast the branch points of N ab (y/s) correspond to those of U ab (y/s). A 
perturbative evaluation of the functions N ab (\/s) is futile since it is at odds 
with the causality constraint, i.e. any finite truncation of the object N ab (\/s) 
will violate the representation (|6]). A summation of an infinite set of terms is 
required to restore causality. This is readily achieved by considering the linear 
integral equation 

roo dw y/s- hm N ac (w) p cd (w) pdbjw) ~ Udb(y/s)] 

^ Wr 7T W - piM W — y/s 



c,d 

There is neither a guarantee that a solution of (pi) exists, nor that a pos- 
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sible solution is unique [42f48] . The question of existence and uniqueness is 
particularly cumbersome in coupled-channel systems. The condition that the 
scattering amplitude may be evaluated in powers of the generalized potential, 
at least close to the matching scale lxm — tun singles out the unique solu- 
tion we are interested in. If the non-linear integral equation (|8l admits such a 



solution it is generated by (12) 



The linear equation (12) leads to a solution always, which is unique for a 
reasonably behaved generalized potential. However, the solution does not nec- 
essarily define a solution of the non-linear equation ([8]). The problem emerges 
when the function detD(-y/s) turns to zero at some energy point on the first 
Rieman sheet. Then the ratio N/D cannot be a solution to pi) since it con- 



tains a pole not present in the original equation. In order to find a solution 



to the non-linear equation the ansatz (11) may be generalized to allow for 
Castillejo-Dalitz-Dyson (CDD) pole structures [48.42J. Physically this points 
to the existence of a genuine bare resonance or to the effect of coupled channels 
not included explicitly. 

Our strategy in the situation where the function det D(y/s) has an unphysical 



zero, is to include one CDD pole in ( 16 ). It turns out that in all cases considered 
in this work, we are able to construct such a solution at least in the energy 
region we are interested in. If the solution found does not satisfy ^ at higher 
energies, one can regard it as a redefinition of the high-energy part of the 
potential through equation 



It remains to specify the structure of the CDD pole terms. We consider the 
more general ansatz [42] 



s - Hm 

y/s - M C DD 

° dw y/s - \i M N ac (w) p cb (w) 

Mthr IT W - LL M W - y/s 



D a b(\/s ) = 5 ab Z. - r ~ lvl R ab 



-E 



(13) 



with a CDD pole mass parameter Mqdd and some coupling matrix R( D \ 
Though the combination (10) and (13) satisfies the unitarity constraint rt7l) 
it does not lead to a solution of the non-linear equation ([81 unless the linear 



equation (12) is adjusted properly. 



The problem arises how to link the CDD pole parameters with the parameters 
of the chiral Lagrangian. In order to establish such a link we consider in a first 
step the case where the generalized potential has a pole below threshold. The 
assumption the pole to sit below threshold is crucial here. This holds in the 
P\\ channel, where the generalized potential has a pole at the physical nucleon 
mass. Note that this is not the case in the P33 channel. Assuming at first the 
absence of the Roper resonance we demonstrate two equivalent ways of solving 
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(Jsj) . Either take the potential in the presence of the s-channel pole and solve 
the original N/D equation (12), or use an effective, pole-subtracted potential, 
and solve a modified N/D equation. An appropriate generalization of (12) is 
derived in the following. We consider the decomposition 



U. 



ab 



g a mg b y/s - hm 
y/s — m m — jj, M 



(14) 



where the generalized potential has an s-channel pole term at y/s = m < /i t h r 
The strength of the pole term is characterized by the coupling constants g a . 



Given the ansatz ( 10 ) together with ( 13 ) the original result based on the equa- 



tions (10, 11, 12) is recovered unambiguously by the following set of equations 



N ab {y/S) = U^{y/S) 



y/s - Hm 



+ E 



with 



dw 



y/s - M C DD 
s ~ m N ac (w) Pcd (w) [Ui( W ) - Ug{yf S )] 



cd J thhr K W - Hm 



W - y/s 



(15) 



R 



{D) _ m - M CDD 



R 



(B) 
ab 



m — p M \ 
Hm ~ Mcdd 



ab 



E 



°° dw m- fi M N ac (w) p cb (w) ~ 



IHbx 71 W — /J, M 



w — m 



(H M - m 



-Yg a mg b 



x ^ f°° dw w- Mcdd N ac (w) p cd {w) 
E / : — ^^9dmg b 



c d , Wthr it w - hm {w-mf 

^(^ M \T r dW N ac{w)pcd{w)Uj{w) 

+ (m- Mcdd) 2^ / 1 w ^ 

^ J Mthr 7i (w- Hm) (w - m) 



(16) 



The merit of the formal result (15, 16) lies in its specification of the CDD pole 



parameters, R^ and R^ B \ in terms of the parameters, g a -,^-, characterizing 



the subthreshold pole term in the generalized potential (14). The CDD pole 



mass parameter, Mcdd-, is irrelevant. By construction, the scattering ampli- 



tude, which results from (10, 13 16), does not depend on the choice of Mcdd- 



We conclude from the rewrite ( T5fi~6 ) that introducing a CDD pole in the N/D 
ansatz may be viewed as a modification of the generalized potential. Subtract 



from the original potential a pole term according to ( 14 ) and then solve the 
system dl3|15|16[ ). 



In a second step we consider the case of a resonance characterized by a pole on 
the second Rieman sheet. In this case the generalized potential does not have a 
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pole. If one assumed a pole at energies above threshold a contradiction with |8j 
arises. Though the decomposition (14) does not exist we can use (16) to relate 
the CDD pole parameters with the resonance parameters. The modified N/D 
equations ( 15|16 ) are well defined for m larger than threshold. The resulting 
scattering amplitude will still be independent on Mcdd if the replacement 
m — > m — ie is used in (16) and if the real part of the integrals in (16) is 
taken. Still the problem remains how to link the resonance parameters with 
the parameters of the chiral Lagrangian, in particular in our case where the 
Lagrangian does not involve any resonance fields. To do so one may imagine 
incorporating a resonance field to the Lagrangian in an intermediate step. 
If the resonance mass is chosen below threshold the identification of U e s is 
defined unambiguously via ( 14 ) . Next we integrate out the resonance field by 
performing a chiral expansion of the resonance pole term as to absorb its effect 
into local counter terms. Finally the case of interest with the resonance mass 
chosen above threshold follows by an analytic continuation in the resonance 
mass. That implies that the pole term in ( 14 ) must be expanded to the chiral 
order considered. This leads to an effective potential of the form 



m — [i-M k=1 \m — H>m J 

as its proper generalization to the presence of a CDD pole. This procedure 
guarantees that the presence of a CDD pole does not renormalize the local 
counter terms to the given chiral order. Though our chiral Lagrangian does 
not involve resonance fields the parameters g a and m have the interpretation 
of a resonance coupling and mass parameter. In our scheme they arise as a 
correlation of an infinite set of higher order counter terms, which we sum in 
response of the causality request. 



2.2 Generalized potential, crossing symmetry and conformal mappings 



At leading order one may try to identify the generalized potential, U a b(-\/s), 
with a partial-wave projected tree-level amplitude. After all the tree-level ex- 
pressions do not show any right-hand unitarity cuts. However, this is an ill- 
defined strategy since it would lead to an unbounded generalized potential 
for which (fS^h does not allow any solution. A remedy would be to restrict the 
integral in Ml) to energies y/s < A s with a suitably chosen cutoff characterizing 
the size of the effective Hilbert space. However, even in this case the gener- 
alized potential needs to be known reliably for energies \fs < A s . A priori 
this restricts the choice of A s to energies where strict converges, i.e. to 
energies where the summation implied by ^ is irrelevant. 

The key observation we exploit in this work is the fact that the solution of 
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the nonlinear integral equation (|8j) requires the knowledge of the generalized 
potential for y/s > /z tnr only. Conformal mapping techniques may be used to 
approximate the generalized potential in that domain efficiently based on the 
knowledge of the generalized potential in a small subthreshold region around 
[ie only, where it may be computed reliably in x?T. To be specific we identify 
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for 
for 
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A function /(-y/s), analytic in a domain Q, can be reconstructed unambigu- 
ously in terms of its derivatives at a point /ig 6 !1, The reconstruction is 
achieved by a suitable analytic and invertible function with £(/i£) = 0, 

which maps Q onto the unit circle in the complex £-plane with |£| < 1. Given 
such a conformal mapping the Taylor expanded function 



OO -L 

Kz) = m-\z)) = Y.^z k 

fc=0 



(19) 



converges inside the unit circle if expanded around z — 0. Thus the series 



fc=0 



(20) 



recovers the original function f(\fs) in its analyticity domain. In contrast to 
a Taylor expansion of the function f(y/s) around yfs = he the convergence 
is no longer limited by the circle touching the nearest branch point. For an 
explicit example and illustration see e.g. 



Given a suitable conformal mapping, £(\/s), we seek to establish a represen- 
tation of the generalized potential of the form 



U(y/s) = U insidc (y/s) + U outsidc (y/s) , (21) 

OO jkTT (C—1(C\\ 

^outsidc(V / S ) = ^ Uk 
k=0 

where we allow for an explicit treatment of cut structures that are inside the 
domain Q. The splitting of the potential into the parts having only poles and 
cuts inside and outside Q is defined up to a polynomial in y/s. Therefore 
C4nside(\/s) can be chosen to approach zero asymptotically. We insist that the 
series for £/ utside(\A) is bounded asymptotically, when truncated at any finite 



t( r \ TT u outsiders \<, ) 

' Uk = Turk 
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Fig. 1. Singularity 
structure of the 
ttN —t- ttN ampli- 
tude (solid lines), and 
domains of convergence 
(dashed lines) charac- 
terizing the expansion 



(21) 



order. Then the expanded potential in (21) is regular enough for the integral 



equation rt8ft to be well defined and amenable to a solution via (12). 



In order to construct suitable conformal mappings it is necessary to discuss 
in some detail the analytic structure of the considered amplitudes. The singu- 
larities of the ttN scattering amplitude are given by the solid lines shown in 
Fig. [I] [5rJfoT)] . Apart from the s-channel unitarity cut extending from thresh- 
old, y/s = rriN + m n , to infinity, there is the u-channel cut defined by the 
intervals 



< y/s < m N — m n , m N - < y/s < Jm 2 N + 2 ml , (22) 

m N v 

the whole imaginary axis and the t-channel cut containing the circle \s\ = 
m 2 N — m 2 , and also the nucleon pole at y/s = imjy. Note that each cut has 
its mirror partner (y/s — > — y/s) because they originate from a Mandelstam 
representation written in terms of s, t, u. 

In contrast to the ttN scattering amplitude its generalized potential is an an- 
alytic function of y/s for energies larger than the s-channel threshold. Thus 
it may be Taylor expanded around the point y/s = \ie- However, the conver- 
gence radius would be limited by the location of the nearest branch point. The 
radius of convergence may be increased if some cut structures can be evalu- 
ated explicitly, such that the nearest branch point of the residual potential is 
more distant. For instance one may consider the u-channel cut implied by the 
nucleon-exchange process explicitly. As a consequence the residual potential 
has its nearest u-channel branch point at y/s = — m^. Similarly the t- 
channel cuts close to the real line may be evaluated in %PT approximatively, 
it would start to build up at order Q 3 . However, there is a natural limit to 
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this construction: the u-channel and t-channel cuts turn non-perturbative at 
some energy necessarily. 

Since we do not know a priori where the u-channel cuts turn non-perturbative 
we insist on that the domain Q, which we seek to specify, excludes the u- 
channel unitarity branch cuts, i.e. 



(— rriN + m n , — m n ) <f_ fl & (— i oo, +i oo) <f_ Q . (23) 

Incidentally such a construction lives in harmony with the crossing symmetry 
of the scattering amplitude. Given the assumption the generalized potential 



cannot be evaluated by means of the representation (21) for y/s < — m 



The decomposition (21) is faithful for energies \fs G Q only. Nevertheless, 
the full scattering amplitude can be reconstructed at yfs < — from 
the knowledge of the generalized potential at \fs > + m n by a cross- 
ing transformation of the solution to ([8]). Such a construction is consistent 
with crossing symmetry if the solution to ^ and its crossing transformed 
form coincide in the region — m n < y/s < + m n . This is the case 
approximatively, if the matching scale fiM in ([8| is constructed properly. For 
mjv — rn n < \lm < tun + ra^ the scattering amplitude remains perturbative in 
the matching interval, at least sufficiently close to yfs ~ hm- With our choice 



fi M = m N , (24) 
this is clearly the case (see also [22J). 

Before providing a first example for a suitable conformal mapping there is 
yet a further issue to be discussed. Since we neglect intermediate states with 
mass larger than + 2 m w one may argue to incorporate a cutoff into ^ at 
\/s < A s ~ ttin + 2m„. This would induce a branch point of the generalized 
potential at \/s = A s . If we insist on the condition 



(A a ,oo)£ n 



(25) 



the effect of higher mass states is part of ?7outside(\/5 ) in (21 ). Thus their effect 



renormalize the counter terms of the effective field theory, as they should. Mod- 
ulo subtle effects from the chiral constraints it is impossible to discriminate 
the effect of a local counter term from the contribution of higher intermediate 
mass states in the generalized potential. 



Though it is possible to use a sharp cutoff A s in (|8j) and construct conformal 
mappings subject to the conditions (23, 25), it is advantageous to use the 
limit A s — > oo in ([8]). This avoids an artificial behavior of the scattering 
amplitude close to the cutoff. Rather than modifying the integral equation it 
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suffices to adjust the generalized potential, U(yfs), such that the influence 
from the region y/s > A s is irrelevant. This is readily achieved by replacing 
^Ca(V^) in (pll) with 



l«A.) 



for 
for 



s < A, 



0. 



(26) 



where the condition £'(A S ) = guarantees a smooth behavior of £a(\/s) at 
the cutoff scale. Though it is not immediate from the nonlinear equation (J8|, 
our claim follows from (12). For a generalized potential that is constant for 
y/s > A s , the contribution of the integral dw in (12) from the region w > A s 
vanishes identically for external energies y/s > A s . 

Consider the dashed line C\ in Fig. [T] which encloses a domain = Qi and 
meets our requirements (23, 25). A conformal mapping of the domain Qi onto 



the unit circle is readily constructed with 



<f)(y/s) +(f)([l E ) ' 

yJ(sA 2 a -(m 2 N -miy) (A^-V) 
s — m 2 N + ml 



(27) 



Note that the construction of (21 ) does not require an explicit representation of 
the inverse function £ -1 (z). The inverse of (27) is not particularly illuminating 
and therefore not given here. The coefficients Uk may be derived in application 
of the chain rule with 



It is not always convenient to map the unit circle onto the largest possible 
analyticity domain of the generalized potential, as we succeeded to do by 



means of (27) in the case of elastic wN scattering. For other reactions finding 



such a transformation can be quite complicated. 



The representation (21) need not converge in parts of the complex plane very 
distant from the physical region. Therefore it suffices to use an universal map- 
ping for ttN — » ttN, syN — > tcN and 7^ — > reactions. 

The transformation function £(\/s) is chosen to be a superposition of the 
function (A s — y/s) 2 and a Mobius transformation, 
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1/2 




A 



2 

ai(A s -y/g) -1 

(a!-2a 2 )(A s - v^) 2 + l 

1 1 

a 2 



(A s - 



2 i 



(As - A c 



Fig. 2. Singularity 
structure of the 
'yN —7- ttN ampli- 
tude (solid lines and 
dots), and domains of 
convergence (dashed 
lines) characterizing the 



expansion (21 ). 



(29) 



For Aq = m\—m\ it maps the interior of the contour Ci in Fig. [Uonto the unit 



circle. Thus it meets our conditions (23 25), in particular it holds £'(A S ) = 0. 



Though not necessary, it is convenient to have available the inverse mapping 
with 



l + x 



\J(l + x) [(1 - x) Oi + 2 a 2 x] 



(30) 



The particular form of the conformal mapping used in this work is irrelevant. 
We constructed various mappings that conform with the physical requirement 
as discussed above. Our final results show a very minor dependence on the 
choice thereof, that can be compensated by a slight change of free parameters. 



We turn to the analytic structure of the 7 AT — > tiN amplitude [HO], which is 
illustrated in Fig. [2j The unitarity s- and w-channel cuts are the same as in 
the case of 7rA^ amplitude, except that the part of the u-channel cut lying on 
the real axis is given by 



< s < 



m N 



m 2 N — m N m n — m\ 



(31) 



vn n + m T 

The t-channel cut from multiple pion exchanges consists of the whole imagi- 
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1/2 

s 




2m^ + m -m . (m 2 - ml) J4m 2 N -m 2 t 

= — * ± i » , 

2 2m t 

m t G (2rn^, 2m N ) . (32) 

In addition there are the t-channel singularities from the one-pion exchange, 
that includes the imaginary axis, the line < s < (mjv — m n ) 2 and a pole at 



m N . We use the conformal mapping (29) with Aq — 901 MeV chosen such 



that the contour C 3 in Fig. [2] touches the edge of the t-channel cut. 

The singularity structure of the 7iV — > 'jN amplitude is presented on Fig. [3] 
The unitarity s- and u-channel cuts are again the same as in the case of 7riV 
amplitude, except that the part of the w-channel cut lying on the real axis is 
given by 



The t-channel cut from the one-pion exchange consists of the whole imaginary 
axis and parts of the circle defined by 



m. 



m N - 



±i 



\JA m 2 N - m 2 , m t G (m.,,2 m N ) 



(34) 



The cut lines of a multiple n-pion exchange include the whole imaginary axis 



and a part of the arcs (34) with m t > nm^. We use the conformal mapping 
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(29) with Ao — 877 MeV, that maps the interior of the contour C4 on Fig. |3| 
onto the unit circle. The particular value for A makes the contour C4 not 
to cross the part of the t-channel cut, corresponding to two and more-pion 
exchanges. 

In the generalized potential based on the tree-level scattering amplitudes (|2j 
111 IBT) only the s- and u-channel nucleon exchange and the t-channel pion- 



exchange contribute to £/i n side(\/s ) i n (21)- The corresponding singularities 



must be separated in 'inside' and 'outside' parts as shown in Appendix |B| 



The 'outside' part of the potential is expanded as £/ utside(\/s ) i n (21)- Given 
the particular choices of the conformal mappings the one-loop expressions 
for the three reactions contribute to U ou t B a e (y/s) only. The order to which 
^outsidc(v / s) must be expanded is discussed in the result sectional The goal is 



to correlate the order of the expansion (21) with the chiral order as to obtain 
a generalized potential accurate in the region /ia/a/s < A s uniformly. Then 
the typical error may be identified in the chiral domain around yfs ~ fiM- 

There is one peculiarity appearing when we consider the one-loop contribu- 
tions to the generalized potential. We use a heavy baryon prescription for loop 
diagrams because it is well established for all reaction channels considered. 
Since the potential must possess only left-hand singularities it is necessary 
to subtract cut structures associated with the s-channel unitarity cut. The 
latter is responsible for the imaginary part of the one-loop expression at en- 
ergies ^/s > m n + mjv- In accordance with ^ we subtract from the one- loop 
expressions the integral 



E 

c,d 



r°° dw 

IfHbxB 7T W 



UW(w)p cd (w)U { 2\w) 



w 
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(35) 



in a given partial wave, where U^(y/s) is the leading order tree-level potential. 



In (35) the reference to angular momentum and parity is omitted. The inte- 



gral (35) is finite since the leading-order potential is asymptotically bounded. 



This is a consequence of the representation (21 ). The one-loop expression after 



subtraction of (35) still has a residual imaginary part at yfs > m n + tjin- It 



diminishes systematically if the order of the truncation in (21 ) and the order of 



chiral expansion is increased. We neglect this imaginary part as higher order 
and decompose the residual loop-expression according to (21). The trunca- 



tion of the sum in (21) is uniform for all chiral moments of the potential, in 



particular for the leading-order term specifying the subtraction integral (35). 
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3 Numerical results 



In this section we present our numerical results for the three reactions tcN — > 
ttN, syN — > ttN and 'yN — > 7 AT, where we assume the ttN channel in an s- or p- 
wave state. We consider the energy region from threshold up to 1.3 GeV, which 
is motivated by the expected reliability of the two-channel approximation 
(the inelasticities in the tcN scattering are small at these energies [51f52] ). 
A detailed comparison with available experimental data is provided. We aim 
for a description of the experimental data set by solving the integral equation 
(|sj) with the generalized potential calculated in perturbation theory based on 
the chiral Lagrangian 0. The generalized potential is analytically continued 
beyond the threshold region and expanded systematically in application of 



suitably constructed conformal mappings as discussed in section 2.2 



3.1 irN elastic scattering 



The irN sector is developed independently from the 7 A channel as it is treated 
to first order in the electric charge, e ~ 0.303. The chiral Lagrangian ([I]) 
provides several free parameters. To order Q the two parameters, / and g^, 
are relevant, where / ~ 92.4 MeV is identified with the pion decay constant 
and gA — 1-26 with the axial- vector coupling constant of the nucleon. Owing to 
the Goldberger-Treiman relation f g^NN = Qa^n one may use alternatively 
the empirical pion-nucleon coupling constant 



^^-13.54, (36) 
4tt v ' 

as input [53J. We do not count the masses of the pion and nucleon, and m N , 
as free parameters. They are assumed to take their empirical values properly 
isospin averaged. At chiral order Q 2 four additional parameters ci, 02,03,04 
arise. There remain the four counter terms proportional to di + d 2 ,d^,d 5 , du — 
di5 that turn relevant at chiral order Q 3 . The parameter d\$ parameterizes 
the degree of violation of the Goldberger-Treiman relation. A study of pion- 



nucleon scattering only does not determine so that we use (36) as input. 
Given the empirical value for g^NN the parameter <i 18 enters no longer as a 
free parameter in our study. 

The generalized potential in (Tsl) is decomposed into an inside and outside part 



based on the conformal mapping (29). While the inside part of the general- 



ized potential is determined by the pion-nucleon coupling constant g^NN, the 
outside part depends on the various counter terms and the order at which 



the expansion in (21) is truncated. As discussed in the Appendix B we per- 
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Fig. 4. Results of the fit for irN S and P-wave phase shifts. The solid curves corre- 
spond to the full Q 3 results, the dashed curves to Q 2 results, and the dotted curves 
to Q 1 calculation. The dash-dotted line in the Pn phase denote a phenomenological 
fit as explained it text. The data are from |51j (circles) and [52] (squares) . 



form a summation of the outside part as is implied by the contribution of the 
u-channel nucleon-exchange cut for a/s > — m^/m^. This controls large 
cancelations amongst the inside and outside parts of the potential that arise 
for large angular momentum. We will determine the truncation order of the 
residual contributions in (21 ) by the number of free counter terms contributing 
to the outside part of the potential. The original potential contains polynomial 
terms in yfs that are unphysical at large energies. In contrast the extrapolated 
potential is bounded at each order in the expansion of (21). Thus, there is an 
issue how many terms in the expansion should be considered. If too many 
terms are included the resulting potential would be unphysically large, even 
though the potential would be bounded asymptotically. Since the outside part 
of the potential is governed by left-hand cuts that are far distant, the expan- 
sion (21) should converge quickly if applied to the full potential. This is nicely 
confirmed by the following phenomenological study. 



In a first step we consider the inside part of the potential as determined by the 
pion-nucleon coupling constant (36) but keep the values Uk in (21) unrelated 
to the parameters of the chiral Lagrangian. We use A s = 1500 MeV, but 
checked that moderate variations with 1300 MeV < A s < 1600 MeV lead to 
almost identical results. Within this phenomenological framework we study 
the relevance of higher order terms in the expansion of the residual outside 
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part of the potential. For the fits we take the energy independent phase-shift 
analysis [52]. The analysis from |5T] is used for comparison. 



The empirical partial waves Sn,Psi and P13 can be well reproduced from 
threshold to \/s < 1300 MeV with only the leading term in the expansion. 
The S31 partial wave requires two terms. The results for these four partial 
waves are shown by solid curves in Fig. [4] and compared with the empirical 
phase shifts [5T1I52] . 



The phase-shifts are obtained as solutions of the linear equation (12). The 
original non-linear equation ^ is satisfied for the Sn and P13 partial waves 
for energies from threshold to \fs < 1300 MeV with a violation smaller than 
1%. The deviation is somewhat larger for the P31 and S31 waves, for which 
we obtain about 10% at the upper end. It is possible to introduce a number 
of CDD poles to get solutions that satisfy the causality constraint more pre- 
cisely, however, this would introduce additional free parameters without clear 
physical meaning. The small deviations obtained we take as a signal of higher 
energy effects not yet controlled in the present work. 

One cannot obtain a satisfactory fit for Pu and P33 partial waves based on 



the linear equation (12) without severe violation of the causality relation (|8j). 
This reflects the presence of low lying resonances in these partial waves, the 
A(1232) and the N(1440) resonances. In the Pu partial wave a numerically 
large contribution of the nucleon pole graph to the potential makes it impos- 
sible for the solution (|8]) to exist without a CDD pole. Similar findings were 
discussed previously in models [54"|f55] also based on a single-channel approx- 
imation. After inclusion of one CDD pole, we can reproduce the P33 partial 



wave with only the zeroth order coefficient in the expansion (21). The Pu 
partial wave is well described in the presence of one CDD pole if the corre- 
sponding potential is expanded to first order. The results are shown in Fig. [4] 
by the solid line (P33) and dash-dotted line (Pu). In both cases the set of 



equations (13, 15, 161) provides solutions of (pi) accurate at the 1% level. 



We return to our effective field theory which relates the expansion coefficients 
in (21) to the parameters of the chiral Lagrangian ([I]). At chiral order Q 



2 



there are four counter terms, Ci, ..,04 which contribute to the two s- and four 
p-wave potentials. There is no contribution to the potentials of higher partial 



waves. The coefficients Uk in (21) posses a chiral decomposition into powers of 



the pion mass modulo logarithm terms. For the p-waves we have to truncate 



the expansion (21) at leading order, k = 0, since any additional term may be 
altered arbitrarily by higher order counter terms. For the s-wave potentials 
the second term with, k = 1, may be considered in addition, since the Q 2 
counter terms determine the slope of the s-wave potential to leading order in 
a chiral expansion. However, it is unclear a priori whether one should go to the 
maximum order introduced above. Since the terms with k > 1 affect the high 
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Table 3 

Parameters obtained from a fit to the empirical s- and p-wave pion-nucleon phase 
shifts. 



energy part of the potential we find it advantageous to consider such terms 
only at an order where they are determined quite precisely. To be on the safe 
side, we restrict ourselves to the zeroth order in the expansion (21) at chiral 
order Q 2 . At the next order Q 3 there are additional counter terms available 
that lead to a precise determination of the slope of the s-wave potentials. 
Therefore it is justified to consider the terms with k = 1 in (21). For the p- 
waves we keep the k = terms only, since the slopes of the p-wave potentials 
receive contributions from unknown Q 4 counter terms. 



We determine the chiral parameters by a fit to our phenomenological poten- 
tials, discussed above. At chiral order Q 3 it can be reproduced identically with 
the exception of the P\\ potential. The numerical values of our parameters are 
collected in Tab.|3j The solid lines in Fig.|4]show the resulting phase shifts. All 
phase shifts coincide with the previously discussed phenomenological phases 
with the exception of the Pu phase, for which some further improvement is 
desirable. An accurate reproduction of the phenomenological Pu potential 



would require the inclusion of the term Uk with k — 1 in (21). Such a term is 
not generated in our approach convincingly. The Pu needs further attention, 
in particular the inclusion of the inelastic ttttN channel (see e.g. p3], [T5] . 

EH]). 



The convergence properties of our approach are illustrated in Fig. [4] by addi- 
tional dashed and dotted lines. The dashed lines correspond to the Q 2 calcu- 
lation. It follows from the Q 3 result by switching off the contribution of the 
one-loop diagrams together with the parameters, di + <i 2 , d 3 , <i 5 , c?i 4 — di 5 . In 



addition the terms with k = 1 in (21) are dropped. The dotted lines show the 
results with c, = 0, in addition, which defines the order Q result. We do not 
change the CDD pole parameters when going from Q 3 to Q 2 or Q 1 results. 
The solid, dashed and dotted lines in Fig. [4] are all obtained with the same 
parameters as given in Tab. [3] 

For all partial waves we obtain a convincing convergence pattern, although a 
few comments must be made. In the Pu partial wave there is a convergence, 
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current work KA86[EE] EM98[57] 


(tyN) 

a [s _] iH 
a< fs+] [ fm l 


0.116 0.130 0.109 ± 0.001 
0.003 -0.012 0.006 ± 0.001 
0.010 0.008 0.016 
-0.058 -0.044 -0.045 


f 7T TV") r Q i 

a [Pn] ] 


-0.082 -0.078 -0.078 ± 0.003 
-0.048 -0.044 -0.043 ± 0.002 
-0.032 -0.030 -0.033 ± 0.003 
0.193 0.214 0.214 ± 0.002 



Table 4 

Pion-nucleon threshold parameters (see e.g. |22j). 

but the Q l and Q 2 results are completely off the data. This is due to a signifi- 
cant cancelation between a large contribution from the nucleon pole term and 
the counter terms. The S31 potential receives a sizeable contribution from the 
first order expansion term, which is absent in the Q 1 and Q 2 lines. 

In order to further scrutinize the consistency of our calculation we provide 
the threshold values of the s- and p-wave amplitudes as they come out of our 
calculation. In Tab. [4] they are compared with the results of different partial- 
wave analyses [5TJ57] . Given the discrepancies among the analyses we find an 
acceptable overall pattern. For the p-wave scattering volumes we agree with 
the analysis of [51] at the 10% level. The s-wave parameters spread most widely 
in the isospin even channel. A direct extraction of both scattering lengths from 
an analysis of pionic hydrogen and pionic deuterium data based on x?T [58] 
gives the following values 



ag_, ; = (0.120 ± 0.003) fm , agjy = (0.002 ± 0.003) fm , (37) 

which are consistent with our results. 

We continue with a discussion of our parameter set as determined from a fit 
to the s- and p-wave pion-nucleon phase shifts and collected in Tab. [3} The 
parameters related to the two CDD poles reflect the presence of the A (1232) 
in the P33 and A^(1440) in the Pi% waves. They allow for an interpretation in 
terms of effective resonance vertices of the form 



off 
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ci [GeV^ 1 ] 


c 2 [GcV- 1 ] 


c 3 [GcV- 1 ] 


c 4 [GeV- 1 ] 


current work 


-1.31 


2.03 


-5.11 


3.59 


K-matrix 


-1.60 


3.25 


-6.34 


3.93 


Q 2 phases [59] 


-0.57 


2.84 


-3.87 


2.89 


Q 3 phases [5] fit 1 


-1.23 


3.28 


-5.94 


3.47 


Q 3 phases [5] fit 2 


-1.42 


3.13 


-5.85 


3.50 


NN phases [60] 






-4.78 


3.96 




d t + d 2 [GeV-' 2 } 


d 3 [GcV- 2 ] 


d 5 [GeV" 2 ] 


(du ~ d 15 ) [GeV" 2 ] 


current work 


0.13 


0.79 


+0.68 


-0.49 


K-matrix 


3.94 


-3.68 


-0.17 


-7.46 


Q 3 phases [5j fit 1 


3.06 


-3.27 


+0.45 


-5.65 


Q 3 phases [5] fit 2 


3.31 


-2.75 


-0.48 


-5.69 



Table 5 



Comparison of low-energy constants in pion-nucleon scattering. 

with the isospin transition matrices normalized by T- Tj = 5ij — Ti Tj/3. 

It is interesting to compare our parameter set with the ones from different 
analyses. A collection of various sets is offered in Tab.[5j where we recall central 
values only. We do not show statistical errors arguing that a comparison of 
different parameter sets is more significant in our case. At chiral order Q 2 a 
strict xPT fit to the ttN phase shift was performed in [22] • The sizes of the 
Ci_4 counter terms agrees roughly with our values, which we recall in Tab. [5] 
for convenience. A more quantitative agreement is expected at higher order 
in the chiral expansion. The parameter sets of [5] are based on a strict %PT 
fit at order Q 3 . The empirical ttN phase shifts at energies y/s < 1100 MeV 
are described. Indeed, the Ci_ 4 parameters agree significantly better with our 
set. In [5] different parameter sets are obtained, the most significant two of 
which are represented in Tab. [5j A further source of information on some low- 
energy constants is a chiral analysis of the nucleon-nucleon scattering process 
[50] . Here only the parameters C3 and C4 are relevant. The values reported in 
[50] differ somewhat from the results of [5]. While the value for C4 is quite 
compatible with the result of [60], the value for C3 is somewhat larger than 
our result. 

The size of the parameter c\ can be related to the pion-nucleon sigma term 
[1]. To the order Q 3 it holds 



a nN = -4 ci ml - 



9 ml 
64 7r ml 



2 

' 9nNN 



N 



+ o(q 1 



(39) 
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Fig. 5. Results of the fit for pion-nucleon s- and p-wave phase shifts using a K-matrix 
ansatz. The parameters are given in Tab. [5j The data are from [51] (circles) and 
(squares) . 



which for our value of c\ would imply a n N = 77 MeV. This value is in contra- 
diction to the recent determination of the sigma term based on unquenched 
but two-flavour QCD lattice simulations [51] . which suggests a much smaller 
sigma term a^N = (44 — 54) MeV. Such values appear compatible with an 
analysis of the subthreshold ttN amplitudes performed to chiral order Q 3 [62J , 
which suggests a significantly smaller value for |ci|. The scattering amplitude 
was reconstructed inside the Mandelstam triangle by means of dispersion rela- 
tions using empirical phase shifts. This may hint at a significant sensitivity on 
how to determine the parameter set and reflect the influence of higher order 
effects. Our work does not shed additional light onto this puzzle. 

A striking discrepancy we observe for the Q 3 counter terms. For instance 
our value for | c?3 1 is about about four to five times as small as the values of 
[5|6] . We traced the source of the discrepancy in the Q 3 counter terms as the 
consequence of rescattering effects. It is not related to the expansion of the 
potential using conformal mappings. If instead of solving the non-linear equa- 
tion ^ we insist on a i^-matrix ansatz as used in [ofrj] the parameter set 
is altered significantly. More specifically, in this case we apply the conformal 
mapping technique to the K-matrix. While the tree-level contributions to the 
K-matrix and the generalized potential are identical, the one-loop contribu- 
tions differ. The K-matrix does not receive a contribution of the second order 



rescattering term (35) nor the CDD pole correction term (17). Only the real 
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part of the one-loop diagrams contribute to the K-matrix. It is decomposed 
according to (21), where we use the same truncation order as in our proper 
approach. The result of a low-energy refit of the counter terms in the K-matrix 
ansatz is included in Tab. [5] by the two rows labeled 'K-matrix'. The values 
are much closer to the xPT fits of [5|6] . Note that the C2 parameter is also 
affected significantly by rescattering effects bringing its value close to the one 
obtained in [5]. The resulting 7rN phase shifts are shown in Fig. [5] Such an 
ansatz, which is at odds with the causality constraint of local quantum field 
theory, is nevertheless able to recover the empirical phase shifts amazingly 
well, with the exception of the P33 phase. The latter is described only much 
below the isobar resonance. 
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3.2 Pion photoproduction 



The generalized potential for the 7iV — > ttN reaction is calculated to chiral 
order Q 3 , where a multipole expansion is performed with the irN channel in 
an s- or p-wave state. At leading order the potential is determined by the mass 
parameters and the pion-nucleon coupling constant (36). At subleading orders 
the anomalous magnetic moments of the proton and neutron, k p and 
probed. They determine the isoscalar and isovector moments, k s and k v , with 



t\j S — tZp ~\~ Hifi ; Kjy — Hhp Kj n ; Kjp — 1.793 . Kj>n — 1.913 ; ^40) 

where we recall also their empirical values. At order Q 3 the four counter term 
combinations d 8 , dg, d 2 o, 2 d 2 i+d 2 2 are relevant in addition. Since we allowed for 
a CDD pole in the Pu and P33 elastic ttN amplitudes, there will be additional 
parameters, g Pu and g[ p 33 , g\ p 33 , characterizing the coupling of the CDD pole 
to the 7./V states. An interpretation in terms of effective resonance vertices is 
provided by 



jpf W P ;( J ) 



£ e s = 7?^ A" lv 75 T 3 N + ^\ A M 75 T 3 (d v N) + h.c 



N* _i_ N* 

6 N* Ks . v T3 N F^ + h.c. 



4m N 2 

f (i) _ m N V3gJ N P Jm N -m A )-2 9 J N P33mA 

J A — ? — 0.14 , 

em A 

W = _ 4V3 S %,m^ ^ 58 
em A 

= -7^ (<-«?;',) -i^ (4i) 



The generalized potential is decomposed into an inside and outside part ac- 
cording to the decomposition (21). The inside part is determined unambigu- 
ously by g w NN- It consists of multiple pole terms at ^/s = m N . The outside 
part contains branch cuts located outside the contour C 3 of Fig. [2] only. It 
is expanded systematically using the conformal mapping (29) with A = 991 
MeV and \xe = itln + m n /2. Like for the elastic ttN potential we find it 
advantageous to perform a particular summation in the outside part of the 
photoproduction potentials. The one-pion exchange contribution provides a 
left-hand cut at y/s < m N — m n that is characterized by multiple poles at 
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d 8 GeV 2 


dg GeV 2 


d 20 GeV 2 


(2d 2 i-d 22 ) GeV 2 




■yn 


7 JV 


7 iV 

ff 2,P 3 3 


3.35 


-0.06 


0.61 


0.05 


-0.33 


0.04 


-0.44 


-0.93 



Table 6 

Parameters obtained from photoproduction data. 



the branch point = 



Though such structures may be expanded 



systematically using (29), the convergence would be unnaturally slow in par- 
ticular at large J. In order to treat the pole structures accurately we consider 
the region — 2m n < yfs <m^ — explicitly in the outside part of the 
potential. For more technical details we refer to Appendix B. In the residual 
outside part of the potential we truncate the expansion (21) at k = for the 
s-waves and p-waves. 



The parameter set is obtained from the empirical photoproduction s- and p- 
wave multipoles. There are in part large discrepancies among different energy- 
dependent analyses [63]fl~8] . Therefore we try to adjust the parameters to the 
energy independent partial- wave analysis from |63j . which is less biased than 
the energy dependent multipole analyses. We fit only the real part of the 
multipoles, since the imaginary parts are given by Watson's theorem |64J. Our 
preferred parameters are given in Tab. |6j 
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Fig. 6. Electric photoproduction multipoles Eq+. The data are from |63J (filled cir- 
cles - energy independent analyses, open circles - energy dependent analysis) and 
[18] (squares) . The solid line denotes the Q 3 calculation, the dashed line - Q 2 calcu- 
lation. 
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Fig. 7. Electric photoproduction multipoles E\ + . The data and line conventions are 
as in Fig. |6j 

We discuss some details of the fit procedure. We find the nN rescattering 
effects important for the real parts of the s-wave multipoles. For the p-wave 
multipoles, which do not have a CDD pole contribution, rescattering effects 
are mostly responsible for generating the correct imaginary part, but do not 
modify the real parts much. The s-wave multipoles shown by solid lines in 
Fig. [6] depend on the particular counter term combination 



d 2 o + 



2 d 2 i — d- 



22 



(42) 



Since none of the other multipoles depend on d 2 o, the parameter combination 



(42) is determined by the empirical s-wave multipoles. Given the spread in 
the different energy dependent analyses a satisfactory description of the s- 
wave multipoles is obtained. 

(3 12) 

The magnetic multipole M 1+ multipole provides a large and dominant con- 
tribution to the cross sections in the A (1232) resonance region [63fT8] . As a 
consequence the empirical error bars for this multipole are very small. The 
multipole depends on the particular parameter combination 



d 8 + 



2 d 2 \ — d- 



22 



(43) 
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Fig. 8. Magnetic photoproduction multipoles Mi +. The data and line conventions 
are as in Fig. [6j 



Since the same parameter combination enters the electric multipole E\+ 



(3/2) 



the parameter combination (43) together with the CDD pole parameters 



9i P33 ' #2 p 33 * s determined by a fit of the E[ 3 _l^ and m[+ 2 ^ multipoles. The 
solid lines in Fig. [7] and Fig. [8] show that a satisfactory description is ob- 
tained. The electric p-wave multipoles p E^_l 2 \ n E^l 2 \ also shown in Fig. 7J 
do not depend on any of the parameters collected in Tab. [6j Nevertheless, the 
empirical multipoles are recovered reasonably by the solid lines. 

There remain two parameter combinations, dg and 2 d 2 i — ^22 ; that cannot be 
determined unambiguously from the energy independent multipole analysis 
[63] . Incidentally, the energy dependent analyses [63UT8] differ most signifi- 

(1/2) (1/2) (3/2) 

cantly in those five magnetic multipoles, P M 1± , n M 1± and M{_ , that 
are left to determine dg, 2c? 2 i — d%i and the CDD pole parameters g P p and 
gp xx . In Tab. |6]we provide our preferred parameter set that took into account 
additional constraints from photoproduction cross sections and Compton scat- 
tering data. The five magnetic multipoles are shown in Fig. [8] and Fig. [9} 

Before scrutinizing in more depth the quality of the multipole description 
we illustrate the convergence properties of our approach. In all Figs. [6]j9] the 
dashed lines show the effect of switching off the contribution from the one-loop 
diagrams and counter terms in the photoproduction potentials. The CDD pole 
parameters are unchanged. For most of the multipoles this affects the results 

(3/2) fl/2) 

by less than 20%. For the M\_ and M 1+ multipoles the difference at the 
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Fig. 9. Magnetic photoproduction multipoles Mi_. The data and line conventions 
are as in Fig. [6j 

highest energy considered constitute a factor of two. Also the p M{_ and 
n Mi_ multipoles prove sensitive to the one-loop and counter term effects. 
This is analogous to our findings for the Pn phase shift, the potential of 
which is a result of subtle cancelations. We conclude that all together there is 
a satisfactory convergence pattern. 

We proceed with a more direct comparison of our results with empirical cross 
section data. This is achieved by taking our s- and p-wave multipoles in the 
calculation of differential and polarization cross sections. Higher partial waves 
are supplemented from two different sources. First we compute the higher 
multipoles within our given scheme but neglect the final state interaction, 
i.e. we identify the generalized potential with the partial-wave production 
amplitude. Note that the one-loop contribution of a strict computation 
does not contribute to any of the higher partial waves for which we neglect 
the final state interaction. Second we take the higher multipoles from the 



energy dependent analysis [18J. In Figs. 10 11 we present differential cross 
section data from several recent and precise measurements [U31l67|l66|ISH| for the 
7P — > ir°p and 7p — > n + n and , yn — > ir~p reactions. Fig. 10 illustrates that our 
multipoles describe the differential cross section for the 7r° production off the 
proton accurately, the effect of higher multipoles being quite small. Recall the 
absence of the t-channel pion exchange process in this reaction. After adding 
the effects of the higher multipoles the same is true for the tt production cross 



sections of Fig. 12 and Fig. 11 The solid lines show the contributions from 
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Fig. 10. Differential cross section for the reaction ~/p — > 7r°p. The solid line corre- 
sponds to our calculation with only s- and p-wave multipoles included. The effect 
of higher partial waves is shown by the dashed and dotted lines as explained in the 
text. The data are from 651. 



the s- and p-wave multipoles. The dashed and dotted lines result upon adding 
up the higher multipoles from our theory and [TS] respectively. The effect of 
using different higher multipoles, dashed versus dotted line, is irrelevant for 
our conclusions. 

We compute the beam and helicity asymmetries of the two reactions 'jp — > ir°p 
and 7p — > tt + n, for which empirical data exist. Their relation to the multipole 



amplitudes is given in Appendix A. 3 The results are shown in Figs. 13] 16} 
After inclusion of higher partial wave contributions all considered observables 
are nicely reproduced in the whole energy region considered. Higher partial 
wave contributions are significant in the production of charged pions. 

Like for irN elastic scattering we scrutinize the near-threshold physics for pion 
photoproduction in some detail. We compute the s- and p-wave threshold 
production parameters. Empirically it is well established that in neutral pion 
photoproduction there is a strong cusp effect at the 7r + n threshold [BUfTU] . To 
obtain accurate results we depart from the isospin formulation and perform 
a coupled-channel computation in the particle basis using physical masses 
for the nucleons and pions. Isospin breaking effects are not considered for the 
generalized potential being estimated to be of minor importance. No additional 
parameters arise. 
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Fig. 13. Beam asymmetry for the reaction 7^ — > ir°p with data taken from [65 1. The 



assembly of the lines is as in Fig. 10 
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Fig. 14. Beam asymmetry for the reaction ^yp —> ir + n with data taken from [68]. 



The assembly of the lines is as in Fig. 10 
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Fig. 15. Helicity asymmetry for the reaction jp — > ir°p with data taken from [66] 



The assembly of the lines is as in Fig. 10 
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Present work 


xPT (Q 4 ) 


Experiment 


E 0+ (vr+n) [l^/m^] 


27.4 


28.2 [9] 


28.06 ±0.27 ±0.45 [72] 


E 0+ (n-p) [l(T 3 /"v + ] 


-31.5 


-32.7 [9] 


-31.5 ±0.8 [73] 


E 0+ (7T°p) [10- 3 /^+] 


-1.12 


-1.16 [73] 


-1.32 ±0.05 ±0.06 [69] 
-1.23 ±0.08 ±0.03 [70] 



Table 7 

Threshold values for the E + multipole in the jp — )■ 7r° p, 7 n — )• 7r~ n and 7 j> — )• 7r° p 
reactions. 

In Tab. [7] we present the s-wave threshold parameters for neutral and charged 
pion production. Our results are compared to the empirical values, which we 
reproduce quite accurately. Tab. [7] recalls also the results of a xPT analysis 
accurate to chiral order Q 4 . While xPT appears well converging for charged 
pion production a less convincing convergence pattern is found for the neutral 
pion production [8]. The order Q 3 value obtained in [8] for the £"0+ threshold 
amplitude is ±0.90 • 10 _3 /m 7r +, in striking conflict with experiment. Only the 
inclusion of the order Q 4 terms lead to a value compatible with the empiri- 
cal data. A precise prediction of the threshold amplitude is difficult in xPT 
since it depends sensitively on Q A counter terms that are not well known. 
The threshold value given in [71] is based on an analysis of the latest MAMI 
data set |70j . It is amusing to observe that we obtain accurate threshold pa- 
rameters already based on the chiral Lagrangian truncated at order Q 3 . We 
reiterate that our threshold values are a consequence of a parameter fit to the 
production data at energies excluding the threshold region. 




CM. 



CM. 



Fig. 17. Near threshold differential cross section for the reaction 7^ — > Tr°p with 
data taken from |75|70] . Shown are results from our coupled-channel theory includ- 
ing isospin breaking effects as are implied by the use of empirical pion and nucleon 
masses. The solid lines correspond to our calculation with only s- and p-wave mul- 
tipoles included. The effect of higher partial waves is shown by the dashed lines. 
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Present work 


xPT (Q 3 ) 


Experiment 


A (7r°p) [W- 3 /ml + ] 


10.2 


9.4 [71] 


9.46 ±0.05 ±0.28 [70] 


P2 (7T°p) [10- 3 K + ] 


-10.7 


-io.o [71] 


-9.5 ±0.09 ±0.28 [70] 


P 3 (tt°p) [10" 3 /m2 + ] 


10.3 


10.6 (ZU 


11.32 ±0.11 ±0.34 [70] 



Table 8 

Threshold values of p-wave multipoles in neutral pion photoproduction. 

Further significant information on the near-threshold region of neutral pion 
photoproduction is available. The s-wave amplitude shows a strong energy de- 
pendence with a prominent cusp structure at the tt + n threshold [70]. Since the 
threshold region shows an intriguing interplay of the s- and p-wave multipoles 



we confront our theory in Fig. 17 with empirical differential cross section of 
the Mainz group directly. Given the fact that we did not fit the parameters to 
those data an excellent description is achieved. 

From the near-threshold differential cross section two combinations of p-wave 
threshold parameters may be extracted. It is customary to introduce the fol- 
lowing three combinations of p-wave multipole amplitudes, 



p cm P\ = 3£i+ + Mi+ - Mi_ , p cm P 2 = 3E 1+ - M 1+ + M X - , 
p cm P 3 = 2M 1+ + M^, (44) 
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Fig. 18. Energy dependence of the p-wave amplitude Pi (l.h.p.) and 
(Pi + P| ± Pi) 1/2 /V3 (r.h.p.) in neutral pion photoproduction. The data are from 
|76j . The solid lines follow from the coupled-channel theory. 
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Fig. 19. Energy dependence of the photon asymmetry in neutral pion photoproduc- 
tion from the coupled-channel theory. The solid lines correspond to our calculation 
with only s- and p-wave multipoles included. The effect of higher partial waves is 
shown by the dashed lines. 

which all vanish at the production threshold with p cm = 0. A complete deter- 
mination of all three p-wave threshold amplitudes requires additional informa- 
tion. For this purpose a measurement of the near-threshold photon asymmetry 
suffices [70J. 

At leading orders in a chiral expansion Pi and P2 do not depend on any of 
the Q 3 counter terms. The threshold behavior of Pi and P2 is predicted in 
terms of the electromagnetic charge, the pion-nucleon coupling constant and 
the masses of the pions and nucleons only [5)74)71] . In Tab. M we recall their 
numerical values from [74] . The chiral corrections of order Q A were studied 
in [71] and shown to be subject to sizeable cancelation effects amongst chiral 
loops and further counter terms dominated by isobar exchange processes. In 
contrast P 3 receives a contribution from the Q 3 counter term combination 
d$ + dg and there is no parameter-free prediction accurate to order Q 3 . 



In Tab. |we confront our p-wave threshold parameters with those of [TOTTI] . 
Though we are in qualitative agreement with the empirical values obtained 
by the Mainz group there is a significant discrepancy. Contrasting conclusions 
were drawn from the computation to order Q 4 in [71], which was able 
to accommodate the threshold values presented in [70]. It is interesting to 
locate the source of the observed discrepancies. For this purpose we provide a 



comparison with near-threshold data of the Saskatoon group [76]. In Fig. 
the empirical energy dependence of Pi is shown against our theoretical result. 
Given the energy dependence of our theory one would expect a threshold value 
of around 10, a value close to our result. This is striking disagreement with 
the value of 9.3 ± 0.09 extracted in [75], based on a different assumption on 
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the energy dependence. In Fig. 18 we present also the empirical constraint on 
the mean p-wave amplitudes. We find an excellent description of the energy 
dependence where the overall magnitude is overestimated by less than 2%. 

We finally turn to the near-threshold photon asymmetry measurement of the 
Mainz group [70]. A direct comparison is not straightforward since an aver- 
age from threshold to 166 MeV laboratory photon energies was performed. 



In Fig. [19] we show the result of our computation at four different energies 
demonstrating a sign change in the asymmetry at energies below the mean 
value of 159.5 MeV in the Mainz experiment. At the mean energy our photon 
asymmetry is about a factor four to five smaller than the averaged value pre- 
sented in [70] . Within our scheme we have no freedom to significantly increase 
that value without destroying the successful description of the photoproduc- 
tion multipoles at higher energies. It is interesting to recall that a negative 
photon asymmetry was obtained also in [77|78] based on a dispersion-relation 
analysis of photo production data. Given such a sign change an average over 
the photon energy depends on the very details of the averaging procedure. 

A significant effect of higher partial wave contributions on the asymmetry is 



illustrated by the dashed lines in Fig. 19 The possible importance of d-wave 



amplitudes in the photon asymmetry was pointed out in [79J recently. 

Since the near-threshold behavior of the photon asymmetry is subject to subtle 
cancelation effects it would be important to extend our analysis to order Q A 
and see whether the observed sign change persists. Also further data taking 
on the photon asymmetry like planned and ongoing at Mainz [80181] is highly 
welcome. 
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3.3 Proton Compton scattering 



We turn to Compton scattering off the proton. Since counter terms start to 
contribute at the order Q A there are no additional free parameters to be de- 
termined in our work. The CDD pole parameters are set already from pion- 
nucleon scattering and pion photoproduction. In the following we present pa- 
rameter free results based on a generalized potential accurate to order Q 3 . 



Since we neglect the intermediate jN states the Compton amplitude reduces 
to the sum of two contributions. There is a direct term and a rescattering 
term, which depends quadratically on the photoproduction amplitude. The 
latter has already been calculated and presented in the previous section. In 
the rescattering part we consider tcN states in s- or p-waves only. Higher 
partial-waves cannot be treated here because those were neither considered in 
elastic pion-nucleon scattering nor in pion photoproduction. We argue that this 
suffices for calculating cross sections in the energy region considered. Higher 
partial waves are largely suppressed by the nN phase space proportional to 
Pcm +1 entering the integral in (|8|. On the other hand for the direct contri- 

3 

_■ 



bution we take into account all J = ~ and J = | waves. The positive and 
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Fig. 20. Differential cross section for Compton scattering off the proton as a function 
of the center-of-mass scattering angle. The data are from [82J (•) and [83J (□). The 
solid lines follow from the our partial- wave amplitudes with J < 3/2. The dashed 
lines show the effect of partial- wave contributions with J > 3/2 as explained in the 
text. 
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Fig. 21. Differential cross section for Compton scattering off the proton as a function 
of the center-of-mass scattering angle. The data are from [83J (□),[68j(A), |84j(o) 
and the lines are as in Fig. [20] 



negative parity partial-wave amplitudes are of equal importance as a conse- 
quence of gauge invariance. We checked that these waves reproduce the Born 
diagrams quite accurately up to energies of about 1300 MeV. According to 
our strategy we have to apply an analytic extrapolation to the generalized 
Compton potential. Using the conformal mapping as detailed in section 2.2 
we truncate this expansion at zeroth order. Terms beyond the zeroth order 
cannot be justified since Q A counter terms may alter those significantly. 

Our results for the differential cross section and beam asymmetry are presented 



in Figs. 20 22 against empirical data. We find agreement with the data taking 
into account some discrepancy of the different data sets. The photon threshold 
region, the pionproduction threshold region, and the isobar region are equally 



well reproduced. This is nicely illustrated by Fig. 23 which shows the energy 
dependence of the cross section at fixed scattering angle 9 ~ 90°. There seems 
to be a systematic undershooting of the backward differential cross section 
for energies \fs < 1150 MeV. We believe, that a natural explanation for this 
effect are missing higher order contributions. Indeed adding to the J < 3/2 
result contribution from additional partial- wave with J > 3/2 as implied by 
the direct term in our scheme the differential cross section is increased towards 



the data in backward direction. This is shown by the dashed lines in Figs. 20 • 



The low-energy physics of Compton scattering is characterized by various 
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Fig. 22. Beam asymmetry for Compton scattering off the proton as a function of 
the center-of-mass scattering angle. The data are from [68] and the lines are as in 
Fig.[20j 



nucleon polarizabilities. For a definition of the polarizabilities see e.g. [85J. It is 
instructive to express the latter in terms of partial-wave amplitudes T± ab (^/s) 
as constructed in Appendix A. 4 and used in (j8|. In our case the indices a and b 
run form one to two, reflecting the two polarizations of the photon. The electric 
and magnetic dipole polarizabilities a p and (3 P probe partial-wave amplitudes 
with J = \ and J = | at threshold, where there are constraints set by crossing 
symmetry. To be specific the following identity holds 



AT ± V* S Wi = m p ) = ° 4 ATZr P 'HVs = m p ) , (45) 
which can be derived from the expressions given in Appendix A. 4 and 



In (45) the 'A' indicates the need to subtract the contributions from the nu- 



cleon exchange processes. In our approach the crossing relation (45) is obeyed 
strictly. This is a consequence of the matching scale hm in ^ being identified 
with the nucleon mass. In turn we recover the one-loop results of [I]: 



a p = — ^ ATZT P 'H™ p ) = ^ 13.0 ■ 10" 4 fm 3 , 

p 16vrm^ +,n V pJ 3847T 2 m^m^ 

ft, - T7T~~ % Ar_y**K) = * l.S-lO-fm 3 , (46) 

16 7rm^ ' 768 n 2 m z N m n 
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Fig. 23. Compton scattering off the proton as a function of energy at scattering 
angle 9 = 90°. The data are from [82J (•), [83J (n),[68](A), [S](o) and the lines are 
as in Fig. 20 The experimental values correspond to scattering angles closest to 



6 = 90° taken in the interval 86° < 6 < 94°. 

which are known to agree well with the experimental values a p = (12.0 ± 
0.6) • 10~ 4 fm 3 , /3 P = (1.9 ± 0.6) • 10" 4 fm 3 [87] (for a recent development see 
[55]). Further polarizabilities probe derivatives of partial- wave amplitudes at 
threshold. Since we expand the generalized potential to the zeroth order in the 
conformal mapping only, we can extract reliably only electric and magnetic 
dipole polarizabilities. 
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4 Summary 



We presented a uniform description of photon- and pion-nucleon scattering 
data up to and beyond the isobar region. The results are based on partial- 
wave amplitudes derived from the chiral Lagrangian formulated with photon, 
pion and nucleon fields. Electromagnetic gauge invariance is kept rigourously. 

Our study is based on partial-wave amplitudes that have the MacDowell rela- 
tions and that are free of kinematical constraints. In the presence of spin the 
derivation of such amplitudes is tedious and presented in this work for pion 
photoproduction and Compton scattering for the first time. The partial-wave 
scattering amplitudes were decomposed into parts with left- and right-hand 
cuts only. Such a separation is electromagnetic-gauge invariant. The part with 
left-hand cuts only defines a generalized potential, which is evaluated relying 
on the chiral Lagrangian truncated at order Q 3 . The part with right-hand 
cuts only is derived as a solution of a non-linear integral equation. The latter 
combines the constraints set by causality and coupled-channel unitary in a 
systematic fashion. The reparametrization invariance of local quantum field 
theory is obeyed strictly. 

The non-linear integral equation does not permit solutions always. The unitar- 
ity constraint implies an asymptotic bound on the generalized potential, which 
is at odds with any polynomial energy dependence as generated by x?T. This 
problem was resolved by performing an analytic continuation constrained by 
the known asymptotic bound. The generalized potential was evaluated in a 
strict chiral expansion to order Q 3 . It followed an analytic extrapolation in 
terms of suitably constructed conformal mappings. The non-linear integral 
equations were solved by N/D techniques, where the presence of two CDD 
poles in the N/D ansatz reflect the physics of the isobar and Roper reso- 
nances. It was demonstrated that the presence of a CDD pole corresponds to 
an infinite set of local counter terms in the chiral Lagrangian. 

Even though the parameter set was adjusted to the scattering data avoiding 
the near threshold region, we obtained results that are compatible with the 
empirical threshold behavior nevertheless. The accurate pion-nucleon scatter- 
ing lengths as measured in pionic-hydrogen and deuterium systems are recov- 
ered. In order to compare with the near threshold data available for neutral 
pion photoproduction we considered isospin breaking effects as implied by the 
empirical pion and nucleon masses. Despite the fact that we reproduce the 
differential cross section accurately we predict s- and p-wave threshold pa- 
rameters that differ somewhat from the values obtained by the Mainz and 
Saskatoon groups. We traced this discrepancy to an additional energy de- 
pendence in the p-wave multipoles that was not considered in the threshold 
analyses of the two groups. As a striking prediction we find a sign change in 
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the photon asymmetry at photon energies below 160 MeV. Further data on the 
photon asymmetry, that do not average over a large threshold region, would 
be highly welcome. Our threshold analysis was completed by a determination 
of the electric and magnetic dipole polarizabilties, for which their empirical 
values are reproduced accurately. 

A detailed comparison of threshold parameters with results from %PT was 
given. In contrast to a strict %PT computation that requires in many cases 
important contributions from the Q A terms, we obtained an accurate descrip- 
tion of threshold observable based on the chiral Lagrangian truncated at order 
Q 3 . The apparent slow convergence of in particular for the neutral pion 
photoproduction process is overcome largely by our resummation approach 
that protects causality and unitarity rigourously. A more convincing conver- 
gence pattern is found in our scheme. 

While in our work we focused on s- and p-waves we provided a comprehensive 
comparison with scattering data up to energies of \fs ~ 1300 MeV. In case 
of pion-nucleon scattering this was achieved by a comparison with the phase 
shifts, which are sufficiently well known from various partial-wave analyses 
and were reproduced quite accurately. In contrast the electric and magnetic 
multipole amplitudes of the pion photoproduction process suffer still from am- 
biguities. A comparison with recent analyses of the MAID and SAID groups 
was offered in this case. Since we observe in part significant discrepancies 
amongst the two groups and our results we provided complementary results 
for differential cross section, beam and helicity asymmetry. We focused on 
the processes which were measured most accurately. For neutral pion photo 
production the influence of additional partial waves was shown to be of minor 
importance. For charged pion photo production the role of higher partial wave 
contributions builds up significantly. A comparison with a representative se- 
lection of high-quality data is offered by combining our results for the s- and 
p-wave multipoles with the higher partial-waves from our theory where we 
neglected the final state interaction. This way an accurate description of the 
data was achieved. The empirical data on Compton scattering off the proton 
were reproduced with the J = \ and J = | partial-waves. The influence of 
higher partial wave contributions was shown to be of minor importance. Dif- 
ferential cross sections and beam asymmetries were reproduced equally well 
from threshold to the isobar region. 
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After submission of our manuscript an improved calculation of aSj and 
a [s+] was announced [89]. The new values are a = (0.122 ± 0.001) fm, 
ag^ = (0.011 ±0.005) fm. 
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A Partial-wave amplitudes and kinematic constraints 



Throughout this appendix we adopt the notation = p^+q^ = Pn+q^ where 
p^ and p^ are the initial and final nucleon 4-momenta respectively, whereas 
g M and g M are the 4-momenta of initial pion (photon) and final pion (photon) 
respectively. The Madelstam variables are s = w 2 , t = (q — q) 2 = (p — p) 2 
and u — {q — p) 2 = (p — q) 2 - In the center of mass frame the energies of the 
nucleons are denoted by E, E, and the energies of initial pion (photon) and 
final pion (photon) by u, Co. For the initial state it holds 



a+r f- m « for ttN f s ~^f+ml f nN 



„2 



N 



for -fN I Vt^ f o r 7^ 



r?2 2 , 2 2 / m ^+Pcm for 

I Pcm f o r 7^ 

with the relative momentum p cm of the center of mass frame. Analogous ex- 
pressions hold for the final state. The center-of-mass scattering angle 9 is 
introduced with 



p» = (E, 0, 0, - Pcm ) , p" = (E, -p cm sin 6, 0, -p cm cos 6) . (A.2) 
A.l Isospin decomposition 



We recall the isospin decomposition for the irN elastic amplitude and the one 
for the photoproduction of the pion (see e.g. [90]). The elastic amplitude T ab 
and the production amplitude T a may be decomposed as follows 



(+) rafe i rp{— ) ■ _abc _ ( 1 _ _ \ rp{\) , f r 1 ^^7-1(2) 



T^.iv = 5" + i e aoc r c = r a r b ) T% + [5 ab - ~ r„ 77, j :/'! v . 

T->a ) • ^aSc 1 /t-i(0) a , ra3 





(A.3) 

where a and 6 are the isospin indices of the final and initial pions respectively. 
The elastic amplitudes with isospin I = ~ and J = | are linear combi- 
nations of X^y • The production amplitudes with isospin I — |, | and 
N = p,n are related to T^ ± ' ) . It holds 
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r^-^to) , 1 i 2 rH _ t (+) , 9T (-) 



s/3 7JV V3 

An "Vdi ^i 7jV ^ 7 iV> 



J 7P - J 7 n _ Y 3 V 77V 1 N J ' 1 tvN — 1 ttN 1 tvN > K 1 ^-^) 

where the amplitudes are introduced with respect to normalized states. This 



is convenient in a coupled-channel framework. The isospin amplitudes of (A. 4) 
are related to the ones in the particle basis as follows 



1 cii 2 

rp rp( 2' I _rp\2> r r > 

1 7T° p— > TT° p g 1 -kN g n N 1 7T° 71— S>7T° 11 1 



2 d"! I cSi 

Ti — _7~ , ^2-'_|__T^^2-' T" 1 



and 



/l (i) /2 (I) _ /l (|) /2 (f) 

J 7p->7r p — y 3 7P y 3 7P ' -^n-s-Tr ™ — y 3 7n y 3 7n ' 



'2 (|) /I (|) /2 (i) /l (l) 

3 V 3 P ' 1 in~^Ti-p — y g J 7« y 3 7 ™ 



A I? 7riV elastic scattering 



The on-shell pion-nucleon scattering amplitude may be decomposed into in- 
variant amplitudes 



T{q,q;w)=F 1 {s,t)+F 2 {s,t)yb 

= F ^){^^} + W){\-^}, (A7) 

where we suppress the reference to the isospin channel and the on-shell nucleon 



Dirac-spinors. The two sets of amplitudes introduced in (A.7) are related to 
each other through 



F± = F 1 ±F 2y /s. (A.8) 
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Since the amplitudes F±, F 2 are free of kinematic singularities, so are the 
amplitudes F± except for the relation 



F-(+yfi,t) = F+(-yfi,t) . (A.9) 

One may adopt the viewpoint that there is in fact only a single invariant 
amplitude F + (y/s,t) which characterizes the scattering amplitude fully. The 
latter amplitude is free of any kinematic constraints. 

The partial-wave amplitudes with definite parity, P, and total angular mo- 
mentum, J, are most economically given in terms of the F± amplitudes. We 
recall from [91122] 



j r L E±m N , 2 nJ-I /-i dcos9 / Pj_i(cosg) 

; _ (E T m N ) 2 F T (^7s,t)\, (A.IO) 



ipL) 



2 



where ± characterizes the helicity of the nucleon. The parity of the amplitude 
t±(\/s) is given by P = ± for J — | odd and P = =p for J — | even. 



The expressions (A.IO) imply the standard normalization of the helicity am- 
plitudes We introduce partial-wave amplitudes T±( that enjoy 
the MacDowell relation 



T J ( +v ^)=T^(-v^), (A.ll) 
and that are free of kinematic singularities and zeros, where we admit con- 



straints at s = as the only exception. Prom (A.IO) it follows that the helic- 
ity amplitudes i±(\/s) suffer from kinematic zeros at threshold and pseudo- 
threshold defined by the conditions 



E ± m N = s + ri f m - ± mN = 0. (A. 12) 

2 v s 

This conclusion is possible only since the invariant functions F±(y/s, t) are 
kinematicly unconstrained. With 



^.) = 2 ^(f)"~\ ± Wi). (A.13) 
t, ± m N \p cm ) 
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we recover the covariant partial-wave amplitudes derived previously in [22] up 
to a factor s J . That factor s J is introduced in the present work as to render 
the phase-space density 



J ( r s ) _ 1 _ Pcm (E±m N ) f pl\ J - 1/2 

P±(VS) - ^Ti(^s) ~ Sirs { s ) ' (A - M) 

constant at asymptotically large sfs. 

We turn to the technical details of the effective field theory. The tree-level 
diagrams (|2]) implied by the Lagrangian density pi) lead to the invariant am- 
plitudes 



F±\Vs,t) = {jj^ (±Vs - m N ) + ^ (q ■ q - (±y/s - m N 



d 3 

+ 772 — 3" {s - m 2 N - m^) (s - m z N - q ■ q) (s - m z N + m 2 n - 2 q ■ q) 

^ J m N 

d x + > 

f 2 m N ^ ' * f 2 m N 



d± + c?2 d. 



+ { - 4 ^ m 2 + ^ 2 (a - - m 2 ) (s - + m 2 - 2 g ■ q) 
v j z j m N 

+ 2 % (q ■ q) + dl t (s-m 2 N -q-q)(q-q- (±y/s - m N ) 2 )X C { + ] 
J 1 r% J 

_ (9a ~2m 2 d ls ) 2 (±y/s - m N ) 2 (J) 

4/2 m „±^ °^ 



+ ^" 2 ^ d ») a ( ^ N {±Vs-m N ) +mN±Vs \ ci n N (A15) 

4 / \ u rrijy J 

where the isospin coefficients C^) are detailed in Tab. [I] The amplitudes 



(A.15) receive contributions from different orders in a strict chiral expansion. 
While the counter terms c... start to contribute at order Q 2 , the counter terms 
d_ contribute at Q 3 . 

At chiral order Q 3 one needs to consider a set of ultraviolet divergent one- 
loop diagrams. The counter term combinations d\ + c^, ^3, ^5 and dn — d\§ 
are required for their renormalization. In the heavy-baryon chiral perturbation 
theory the corresponding loop diagrams were computed in [5] . The results were 
expressed in terms of the frame-dependent amplitudes, g and h, that arise 
naturally in the center-of-mass frame with two-component nucleon spinors, 



_ E + m N 

J- cm — 



(g + i a ■ (q n>ont x h j . (A. 16) 



51 



with the Pauli matrices erj. The invariant amplitudes F± of (A. 7) follow with 



F + = g — hp 2 cm cos 9 , F_ = — h (E + m^) 2 ■ 

We recall from [5] the explicit result 



(A. 17) 
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+ \7^( 2m ' - *X rf i + d *) + -p" rf 3 + — p^ 

-log 
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-t 



si 
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24irf i u 2 



(A. 18) 



where the (i's are the scale independent renormalized coupling constants ac- 
cording to [5]. 



At order Q 3 it is legitimate to add up the loop contribution (A. 18) to the 



tree-level contribution (A. 15) via the relation (A. 17) provided the substitution 
rules 



g A - 2 m 2 dig 



gL ->■ 



m N 



(A.19) 
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are used in (A.15). In the one- loop result (A. 18) we replace qa — > f g-wNN /^n 
for convenience. The loop contributions are specified in a manner such that 
it is justified to use the empirical pion-decay constant with / = /„. in the 
tree-level expressions (A.15). 
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A. 3 Pion photoproduction amplitudes 



We decompose the on-shell pion photoproduction amplitude into four invari- 
ant amplitudes, where we will suppress reference to the channel specifics (see 



(A.3)). There are different choices possible. Chew et al. constructed in [90] the 
CGLN amplitudes, A(s, t),B(s, t), C(s, t) and D(s, t) that are free of kinematic 
constraints [38.92J. Since we are interested in a partial- wave decomposition of 



the invariant amplitudes it is advantageous to use a different set showing the 
MacDowell symmetry [36] in addition. Like for elastic pion-nucleon scattering 
it is possible to introduce invariant amplitudes F^(y/s,t) and F^(yfs,t) with 



Fr(+^S,t)=F t + (- v r S}t ) 



(A.20) 



that are free of kinematic constraints. We recall the definition of the CGLN 
amplitudes and relate them to the F^_ 2 (y/s,t) amplitudes 



i T^q, q; w) = 75 7 M 4 M s , t) ~ 2 7s [ (q • q) P» ~ (P ■ q) %\ B(s, t) 

+ 75 [{q-q)iti- iqj\ C(s,t) 

+ 2 lb [{P ■q) 1 ^-iP fl -m Nl ^i] D(s,t) 

= j:{Ft(V~s,t) L$>+> (q, q,w) + Fr(^, t ) (q, q, w) } 

L£> ± \q,q,w)= l5 ( 1 -± 



1> 



2 IJl 



7/,! 



L^\q,q,w)= 1 J 1 -± iV 



2 2Js 



{q ■ q) 7^ - i q^j , 



(A.21) 



where P M = (p M +p fi )/2. The two representations in (A.21) are equivalent in 
the presence of an on-shell initial nucleon spinor, p 2 = m 2 N and 



Ff = A - B (q ■ q) - D (m N ± yfi) , 



F 2 ± = C - D - B (m N ±y/s). 



(A.22) 



^From (A.22) it follows our claim that the amplitudes F\ (y/s , t), and F 2 (y/s , t) 



are free of kinematic constraints but (A.20). 



The partial-wave amplitudes with definite parity, P, and total angular mo- 
mentum, J, corresponding to anti-aligned and aligned helicities of the initial 
nucleon and photon are readily expressed in terms of the invariant amplitudes 



F 12 introduced in (A.21). We derive 
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(A.23) 



except for J = | with £J. 2 = 0. The parity of the amplitudes i±i(\/s) an d 
^ 2 (\/^) i s given by P = ± for J — | odd and P = =F for J — | even. The 
electric and magnetic multipoles [§Uf§3] can be expressed in terms of helicity 
partial- wave amplitudes t^^y/s) and t± 2 (\/~s) as follows 
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where we made explicit the isospin structure. Note that the conventional elec- 
tric and magnetic multipole amplitudes are introduced with respect to isospin 
states that are not normalized. Since the partial-wave helicity amplitudes are 
defined with respect to normalized states the factors and V| arise. 

The differential cross section, da/dfl, the beam asymmetry, E, and the helicity 
asymmetry, d (03/2 — G1/2)/ dSl, are expressed conveniently in terms of helicity 
matrix elements (see e.g. [SI])- It holds 



da p c 



•'cm 
'cm 



E dn - pV {HspHsA ~ HnHd) 



'cm 



dn 2 p. 

e — = i 

dn % 

d<y 3 /2 _ dai/2 _ P. 
dfl dQ p, 

with x = cos 9 ancS 



(|-^iv| 2 + \Hsa\ 2 + \Hsp\ 2 + \Hd\ 



^ (\H SP f + \H D f - \H N f - \H SA \ 2 ) , (A.25) 

'cm 



Hn = ^ cos 9 -Z - tQ (P' J+ ,(x) - P'^x)) 



H sa = sin °- J2 (tii + tl tl ) (P' J+ M + P'jM) , 



rriN — sin # cos ^ 



f d = ^ T i!^fL^(4 2 + ^)(p; +1 (x) + p; 1 (x)). 



.1 



We construct partial-wave amplitudes, T± A (y/s) and 7± 2 (a/s), that are free 



of kinematic constraints. According to (A. 23) the helicity partial-wave ampli- 



tudes have zeros at thresholds and pseudothresholds 



s + m 2 N - ml 



E±m N = — ^= 7L ±m N = 0. (A.26) 

2 V s 



Taking into account (A.l) and 



PcmPcm COS 6 = - + EE - 171% , (A.27) 

2 Eq. 3.1 of [94J misses a factor I in the expression for Hsa in front of the En+i^ 
term. 



56 



it is readily seen that the two amplitudes 




2m N t±,i(Vs) 



E±m N (p cmPcm ) 
2m N t J ±2 (^fs) 



E±m N (p c 



(A.28) 



are free of kinematic constraints but one exception. From (A.23) it follows 
that the amplitudes t± i(y/s) and t± 2 (\fs ) are linear dependent at u = p cra = 
even after pulling out the common phase-space factor. In that limit the 
two partial-wave amplitudes are characterized by the terms proportional to 
Pcm Pcm COS 9F 2 T (^,t) in (|A.23[| 31 



It follows 



m N 1 2 J-l ±J l f-\ 
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(A.29) 
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Owing to (A.29 ) it is straightforward to derive the following set of partial- wave 
amplitudes 
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ba 



which are free of any constraints. As in the case of vrA^ scattering the factor of 
s J was chosen so that the jN phase-space density is asymptotically constant 
(see (191)). The MacDowell relations 



TtA+ 



T1J+ 



+T1J- 



(A.31) 



hold. 



We turn to the details of the effective field theory. The Lagrangian density ([I]) 
leads with (E A.21) to the tree-level expressions 



3 The amplitudes F±(-\fs, t) and F^(y/s , t) possess dynamical singularities at the 
7 AT threshold: while the s- and u-channel nucleon exchange processes are associated 
with pole terms of the form l/(s — m 2 N ) and l/(u — rn%), the pion-exchange process 
leads to a term l/(i — m^). These dynamical singularities imply a singular behavior 
of the partial-wave amplitude at u = p cm = 0. They appear, however, only in the 
pole graphs which have to be treated explicitly. 
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with the channel dependent coefficients C and C given in Tab. [2] The tree- 



level expressions (A.32) receive contributions at different orders in the chiral 
expansion. 

At chiral order Q 3 a set of ultraviolet divergent one-loop diagrams contribute. 
In the heavy-baryon chiral perturbation theory the corresponding loop dia- 
grams were computed in [8"|fTU] . Their renormalization requires the presence 
of the Q 3 counter terms ^8-9,20-22- The results of PITD] for the •yp — > n° p 
production were expressed in terms of the electric and magnetic multipoles 



(A.25). The one-loop expressions for the yp — > tt + n reaction were presented 
in [10] in terms of amplitudes, that arise naturally in the center-of-mass frame 
with two-component nucleon spinors [HD] • We introduce such amplitudes in a 
notation tailored for our purpose, 



E + m N E + 7777V 

\^r^mV 



aft -i{a-q„) (a x q y ) f 2 



+ (<? ■ <? 7 ) Qn /a + (<? • Qn) q-n , (A.33) 

where we use a different overall normalization and pion and photon three 
momenta that are not normalized. The Coulomb gauge with vanishing zero 



component of the photon wave function is assumed in (A.33). The original 
amplitudes of Chew, Goldberger and Low 



[90] possess kinematical zeros at 
vanishing photon or pion momentum. In part such kinematical constraints are 
a consequence of using normalized photon and pion momenta in the definition 



58 



of the amplitudes. The amplitudes f 2 and fs will be shown to be free of 
kinematical zeros at q n = or = 0. This is not the case for the remaining 
amplitudes fi and f± as can be seen by expressing the invariant amplitudes of 



(A. 21) in terms of the center-of-mass amplitudes 
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Since the invariant amplitudes F 1 and F 2 were shown to be free of kinematical 
constraints, it follows that the amplitudes f\ and must vanish at \/s = m^. 

^From PflQ] we reconstruct the four center-of-mass amplitudes 
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2vr 



CO 



2 + in 



Pen 
id 



C- 



(A.35) 



with renormalized and scale-independent counter terms cf in the convention of 
[8]. The channel dependent coefficients C... are given in Tab.[2j The kinematical 
constraint on the amplitude f\ at y's = m N is incorporated by pulling out the 



overall factor u 



[s — m 



N 



)/(2 



At order Q 3 it is justified to add up the loop contribution (A.35) to the tree- 
level contribution 



(A.32 



via the relation (A. 34) provided the substitution 



rules (A. 19) are used in (jA.32|). In the one-loop result ( |A.35[ ) we replace qa — > 
f 9ttNN /m>N- 
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A. 4 Photon-nucleon scattering 



We decompose the Compton scattering amplitude into invariant amplitudes. 
Bardeen and Tung (92] constructed six amplitudes Af^ 6 (s, t) that were proven 
to be free of kinematic constraints 



Vfa q;w) = lb Af T (s, t) Lf , Q, = , P, = , 

i=l ~ 

L^ = Q 2 g^-2Q^Q", 

Lf = | g 2 (rfQ rf - rfQ rf) - (P ■ Q) (Q^ -f + Q" Y) 

+ ®{Q»P V + Q u Pn, 
Lf = (m N $ - P ■ Q) g» v + i Q 2 o» v - % a au Q a Q»-i Q a Q u 

- m N (Q» Y + Q u Y) + (<5 M P v + Q" P M ) , 

tit = Q 2 (Y P v + Y P 11 ) -(P-Q) Y + Q u Y) 

- $ (Q M P u + Q v F") + ((P • Q)<$ - m N Q 2 ) g» v + 2 m N Q M Q u , 
Lf = Q 2 P^ P u -(P- Q) {Q» P u + Q u P") 

- | (P 2 Q 2 - (P • Q) 2 ) <T + P 2 Q» Q v , 

tit = -\{ p - Q) (Y p u + Y P") + l(P-Q) (Y$ Y-YQ Y) 

- | m N Q 2 a» v + | (m N (P ■ Q) -QP 2 ) gT +Q (Q M Q" + P M P") 
+ | (Q" 7" + Q" 7 M ) - I ™jv (<5 M P" + Q" P M ) 

+ |m JV (7 ai/ g a g^ + im i v(7^g a g^. (A.36) 

For our purpose it convenient to introduce an alternative set of invariant am- 
plitudes, F^_ 3 (y/s,t), for which the MacDowell symmetry |36j is manifest and 
therefore more transparent expressions for the partial-wave helicity amplitudes 
arise. We construct amplitudes which enjoy the relations 



F r(+V~ s ,t)=F t + (-V~s,t), (A.37) 
and are free of kinematic constraints. We find that the following decomposition 
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+ (([w-q) g ua - w u q^j {q a l^ 



(A.38) 



meets our requirements. For on-shell conditions the amplitudes F^L 3 (y/s, t) 
are linear combinations of the amplitudes of Bardeen and Tung [92] . It holds 



if 



m N ± 2 yfsrriN 



i m 



BT 



(m N =F s/sf 



A 



BT 



N 



1 A BT T Vl A BT 



4 



6 J 



iV 



± ± 2 v/s miv + 2 (g ■ g) - s BT 

b 2 = A l 



N 



m 2 N =f 2v/s mAr - 2 (g • g) + s BT 1 BT 



i rriN 



A ai — - A 1 

A 2 ^ A 3 



=F 



.4 



+ 



??? 



A? 



s m 2 N =f 2 s m^v + s 2 — 2 (g • g) s 



4 BT 
A 6 > 



iV 



2m N 1 2 2 



mN^2y/S aBT 

Am N 5 



(A.39) 



We observe that the transformation coefficients, expressing the Bardeen and 
Tung amplitude in terms of the Fi_ 3 amplitudes, are finite everywhere but 
at yfs = 0. Furthermore the determinant of the transformation (A.39) equals 
(y/ss)/ 128. It follows our claim that the amplitudes, F^L 3 (y/s,t), are free of 



kinematic constraints but (A.37). 



The helicity partial- wave amplitudes with definite parity, P, and total angular 
momentum, J, are conveniently derived from the invariant amplitudes intro- 



duced in (A.38). We use a convention consistent with the photoproduction 
amplitudes 



±u 7 



_1 dx 



2 EF? + 2uj Ff ± (2 s - m% ± m N y/s ) F, 



2m N J-i 2 

E F^ + 2co ± m N (m N =f \fs ) F 3 T 



+ 



-2£F 1 ± -2wF 1 T + 4 v / si ? 2 ± T (mfv ±m N ^7s-2s) F 3 T 

+ x ( - 2 E F^ - 2 u Ff ± m N (m N ± y/H) F, 
±co 2 dx 



2 J- 
2 J+3 



Am N J-i 2 
+ x ( 4 mN F i ~{rn N T^f F 3 T 



\( K -Am N F? + (m N ±yTs f F 3 d 
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+ 



+ x VIS ( 4 F i m * - ^ ( m N ± v^) 2 ) 

t J ±j22 (Vs ) = ^J^y{ [tS ( 2 <*> F ? ~ 2 E F ? + fa* ± 
&(2EFf-2uF?-yfi(m N ^y/s)Ff) Pj_i 



I ™ 2 J +3 
* r X 2 J- 



- x ¥jVz ( 2EF ? -2uF^ - y/~s(m N ±^~s 



P J+h (x)\, (A.40) 



with x = cos 6 and P = ± for J — | odd and P = =F for J — | even. Following 
[95] the differential cross section and photon asymmetry are expressed in terms 
of helicity matrix element^"] 



da 

dn 

da 

dn 



2 | 2 + 2|03| 2 + 2|0 4 r + |05| 2 



^(I0i| 2 + I0 

ft ((01 + 5 )0^ + (02-06)0 



+ I06| 2 ) 



4; 5 



with 



(A.41) 



m N 
m N 



An y/s 



m N 



cos ^ E «n + <n) {p'j+dcose) - p'j_ h (co S e)} 



sin 9 E «n -<u) {^ + i(cos^) +P;_i(cos 



sin # sin 



t - 12 t +,12 



{p; + i(cos^) + p;_i(cos^)} 



sin # cos 



4vrv^ x /(J-i)(J + 



m N 



2 cos 



3 e 



3 

2> J 



r E «i2 + <i 2 ) {^ + i(cos^) - p;_i(cos#)} 

E (*-,22 + 4, 22 ) { " 3 (J " §) P^COSfl) 



4tT^ (J-I)(J + 

+ ( j - 1) p; +1 (cos 0) + 2 (p;% (cos 0) - p;: 2 (cos 0)) } 



2sin d 



47Tv^ (J-i)(J + 



3 

2^ j 



E(<22-<2 2 ) {3 U " §) PU (COS, 



We utilize a different phase convention as compared to [95) . The expressions for 
differential cross section and photon beam asymmetry are nevertheless the same. 
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COS I 



+ 2 (P"J_ ! (cos 9) + P"J 3 (cos 



(A.42) 



In order to obtain partial-wave amplitudes that are free of kinematic con- 



straints we apply the transformations derived in (A. 13 A. 30) to the helicity 



partial-wave amplitudes of (A. 40). We introduce 



T lab(Vs) = 2m N y 



s 

Pcm , 



(A.43) 



E 

c,d 



Pcm 





2 J— 1 rrtj 
2J+3 pen 

1 



*i«i(Vi) 



2J+3 p c 



which implies the photon-nucleon phase-space distribution 



db 



p! n ^/S) 



-a 



8vr ^ 



p 2 



J ~ 5 m^VPcm\ 



J+f 



2. "IjVPcm 
S 



i+ J+| , 



/ 



(A.44) 



Inserting (A. 40) into (A.43) we confirm by explicit calculations that indeed the 



partial- wave amplitudes 7± ab (y/s ) are kinematically unconstrained. Moreover, 
the MacDowell relations 



for a 



-Ti^-yfi) for 



(A.45) 



are satisfied. 



We turn to the specifics of the effective field theory. For the photon-proton 
amplitude the Lagrangian density Q and (|5| imply the tree-level expressions 



, ± = e 2 (g A -2 m 2 d is) (m N T /s) 2 e 2 /s 

1 (4tt f) 2 m N (ml - t) m N (m 2 N - u) (m N ± y/s ) 
e 2 K p jm 2 N jm N ± y/s) + s jm N T y/s)) e 2 /t 2 (m N ± y/s) 

m 3 N (m 2 N -u) (m N ± y/s) Am 2 N (m 2 N -u) 
,± = e2 (gA ~2m 2 rf 18 ) ((mjy =F V^) 2 + 

2 ~ (47r/) 2 m^(m2 -f) 

e 2 (mjy (mfy - s) ± -y/i (mfy - tt)) 
m N (m 2 N - u) (m N =F y/s) [m N ± v^) 2 
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+ 



+ 



e 2 ((k p + 2) m N ± K p y/s) 2 
4mfy (m N =f a/s )(mjv ± y/s ) 2 
e 2 k p (2 m A N + s mfy =F 2 \/i sm^ — us) 



m% (m 2 N - «) (mfy - s) 



,± e 2 (c/ A - 2 m 2 d 



18; 



+ 



(27r/)2(m 2 -t) 



+ 



4e 2 



(m^ - «) (m 2 N - s) 



2 e 2 k p (m 2 N =F y/srriN + 2 s) 



which contribute at different orders in a chiral expansion. 



u lmi T - s 



(A.46) 



A systematic computation at order Q 3 requires the evaluation of one-loop 
diagrams. The relevant loops were established by Bernard, MeiBner and Kaiser 
PQ in terms of six center-of-mass amplitudes 



TZ = E + mN \ 5 ij Ai + q i q J A 2 + % e ijk a k A 3 + i e Hm a k q l q m 5 lJ A 4 

Z 777 at 



+ i °k [e«z q q 3 ~ e kji q l q*]A 5 + i a k [e kU q L q ] - e kjl q l q'} A e j ,(A.47) 
within a two-component nucleon spinor formulation. The Coulomb gauge is 



assumed in (A.47) with vanishing zero components of the photon wave func- 



tions. The invariant amplitudes of (A.38) are unambiguously reconstructed 



from their center-of-mass amplitudes, Ai_ e (u,t), with 



a l2 



1 / nix + y/s x " 



s \ 2 m at 



\ ^ 6 

J j=l 



/* + m N -y/s 

a -i o — fli i — ~. , — . „ . d 



(m N + oj (x - 1)) 



l! 11 ~ (m N + yfs) 2 ' " 15 " m N + y/S 

m 2 N ( m N — \/s) (x — 1) — mjy s (x — 7) + y/s s (x — 3) 
4y/s(m N + y/s) 
+ _ (x - 1) + y/s (y/s - 2 m N ) (x + 1) + 

14 — a 7= 5 a 16 



4^ 




771 jy + U) 






= a; (x — 1) 


(mTv — 2 a/s ) (x — 1) 


— s (x — 3) 



-(mjv - y/s) , 



a 13 



a 12 ~ A rz i a 16 



2 (m N + u) y/s 



4^ -ib 

N ' ' 



mil (x — 1) — -v/i (y/s + 2mjy) (1 + x 



1-14 



fl 23 ~~ a 21 



4v/i 



777 TV S (x — 3) + 77% (x — 1) + S2 (1 + x) — m N y/s (1 + 3 x) 
Au) s (mjy + 
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+ m 5 N (x — l) 2 + m% y/s (1 + 2x — 3x 2 ) + 4m 2 N y/s s (x — 1) x 
a 22 ~ 8y/ss(m N + y/s) 

^ y/ss 2 (3 + 2x-x 2 ) -m N s 2 (13 - 2 x + x 2 ) 
8 y/ss (m N + y/s ) 
+ m% (x — l) 2 + Am N y/s s (1 + x) — s 2 (1 + x) 2 

a 24 = 



'S s 

-4 m 3 N y/s (x 2 - 1) + 4 m 2 ^ s (1 + x 2 ) 

____ ? 

8 m% s + m% (x - l) 2 - 3 m% y/s (x 2 - 1) 



a + = 

25 Ay/ss(m N + y/s) 

-y/ss 2 (x 2 - 1) +4m 2 N y/ss (2 + x 2 ) -m N s 2 (5 - 2x + x 2 ) 
4 y/s s (m N + y/s ) 
s (x — 3) + (x — 1) + y/s s (1 + x) — m 2 N y/s (1 + 3x) 

2^ 

8 m N s E + m% (x - 1) + s 2 (1 + x) - 2 m 2 N s (2 + x) 
21 ~ 8u 2 y/s~s ' 

4m N y/s (m 2 N + 2 s) (x - 1) + (x - l) 2 



a 26 



a 22 



°23 



Oo/i — 



8v^s 

s 2 (-3 - 2 x + x 2 ) - 2 s (-3 + 2 x + x 2 ) 

8 77T- 7v s -E + (x — 1) + s 2 (1 + x) — 2 s (6 + x) 

8 uj 2 y/s s 

m% (x - l) 2 + Am N y/s s (1 + x) + s 2 (1 + x) 2 



8yfss 

^ -2m 2 N s(-l +4x + x 2 ) 



8 y/ss 

_ 2m 3 N y/s(x - 1) +m% (x - l) 2 + 2m N y/ss(3x - 1) 

125 ~ 

s 2 (x 2 - 1) -2m 2 N s(x 2 + 5x -2) 

+ 47^ ' 

8 m N s E + m% (x - 1) + s 2 (1 + x) - 2 m 2 N s (2 + x) 

2 S (jTljy - y/s) 

+ + _ 2mjy + _ (x - 1) - y/s m N (x - 3) 

° 33 " a31 "c(m^ + v ^) 2 ' a32 " ^(m w + v ^) 

+ _2m N + m 2 N (x - 1) - y/sm N (x + 1) 

°36 ~~ ~ > a 34 



°26 ~~ 



a 35 = 



u ' \/s(m N + y/s) 

2rriN ijn N — y/s) (x — 1) 



y/s (m N + y/s) 
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m N 



m 



N 



[X 



1) + y/sm N (x - 3) 



A 33 



A 36 



'35 



'31 



UJ- 



'-32 



m 2 N (x — 1 



A 34 



4mjv 
- m N ' 

2 m|r v^s — m 3 N (x — 1) + s m N (x — 1) 

UJ s 



s (m N + y/s) 
sniN (x + 1) 
s(m N - v's) 



(A.48) 



with x = cos 9 and uj,E as in (A.l). We specify the one-loop contribution to 
the invariant amplitudes F^_ 3 in terms of (A.48) and the one-loop expressions 



for the center-of-mass Ai-q(u, t) amplitudes. From |1] we recall 



W = yjm% -u 2 z 2 + t(l- z) 2 x(x-l), R = \jml + t(l- z) 2 x(x- 1) 
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(A.49) 



At order Q it is justified to add up the loop contribution (A.49) to the tree- 
level contribution (A.46) via the relation (A.48) provided the substitution 



rules (A. 19) are used in ( |A.46[ ). In the one-loop result ( |A.49[ ) we replace qa — > 

f 9ttNN /m>N- 
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B Separating singularities of tree-level diagrams 



In this appendix we show how to derive the representation (21) for the nu- 
cleon s- and -u-channel and the one-pion-exchange diagrams. The treatment 
of the nucleon s-channel poles is straightforward. While the pole term at 
y/s = m N is part of £A ns ide(\/s ), any residual background contribution is part 
of U ontsidc (y/s). 

Consider the w-channel nucleon exchange in irN elastic scattering. It takes the 
generic form 



dcos 6 X(\/s, cos 9) 

x ~^T u-m% ' (R1) 

where the function X(y/s, cos 6) reflects the specifics of the partial- wave consid- 
ered. It is a polynomial in cos 9. The contribution to the generalized potential 
is derived in application of the following identity 



dcos 6 1 /* — 1 ds' f — 1 ds' 

+ 



2 u-m 2 N J 4p2 s > _ s J 4p2 (V) s ' - s ' 

1 -oo (A")2 

ml T - ml 



A+ = ^m% + 2ml, A N = ^- , 

niN 

where we consider an s-wave type angular average with X(^/s, cos 6*) = 1 for 



simplicity. A translation of (B.2) from the variable s to yfs is performed by 
means of the identities 



' A ' ^ (a(s') + 0(s') ^~s ) ds' + } N (a(w 2 ) + (3(w 2 ) w) dw 



W — \/S 



N ' ' 2 



a(w 2 ) + (3{w 2 ) w) dw 
w - Js 



(B.3) 



where a(s) and f3(s) are analytic in the integration region. Any contribution 
of the w-channel nucleon exchange to the partial-wave potential can therefore 
be represented as, 
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dcosO A(a/s, cos#) / 7(u7,x ) 



y 7(w,x e J dw . r 

I ~A 9~~? 9\ r ' M outside(V S J 

j Av 2 [w 2 ) w — \/s 



2 it — m 2 N J Ap% m (w 2 ) w 

A N 

2ml w 2 - w 4 + (m 2 N - mp 2 

Xe ~ A 2 2 / 2\ ' ^ A 

Aw 2 p 2 cm (w 2 ) 

where the function u outS idc(\/s ) has no singularities inside the domain enclosed 
by the contour C 2 of Fig. [T] The value xg corresponds to the specific scattering 
angle determined by the condition u = m 2 N . The cut from < \fs < A]v in 



(B.4) contributes to the inside and outside part of the generalized potential 



as seen in Fig. [T] Nevertheless we include both parts explicitly as implied by 



(B.4). This amounts to a well defined summation of an infinite set of terms 
in t/ ou tsidc("\A)- F° r higher partial- waves there would be a large cancelation 
of the inside and outside parts, which we avoid by the suggested summation. 
Such a procedure is justified, since there are no analogous cuts stemming from 
higher order contributions. 

We turn to the — > nN reaction. In this case we need to consider the 
nucleon u-channel and t-channel one-pion exchange. It holds 



dcos9 1 r — 1 ds' 



u-m 2 N J 4p cm (s')p c 

—00 

arctan ( 2 m^ 



\Pc m (m 2 N )\ 

2m 2 N -ml J TUN 
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\Pcm(m%)\ s-mj 



dcosO 1 f — 1 ds' 



ml J 4p cm (s')p c 

—00 

(mjv-m^-) 2 

+ 



ds' 







4j9cm(s / )Pcm(s / ) S 1 



arctan ( 2 m N b " cm ^ }l 



m N 



\Pcm(m 2 N )\ s-m 2 N ' 



^2 



(B.5) 



s - (m N + m n ) 2 2 (s — m 2 N ) 2 



{ m N -m n} ) — , rtaBV . y As 



From the various contributions in (B.5) only the pole terms, singular at yfs = 
771 jy, will be part of inside )• Like for the case of elastic ttN scattering 
it is advantageous to perform a summation in {7 utside("\/s ), as implied by 
the cut contribution from m^ — 2m n < yfs < m^ — m n in the one-pion 
exchange contribution. This term leads to numerically large pole structures at 
y/s = m N — m n . The lower bound at y/s = m N — 2 m n is somewhat arbitrary, 
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but suffices to identify the important pole terms. We assure that our results 
do not depend on that specific choice. As in the case of the nucleon exchange 
in 7iN elastic scattering no such cut structures appear at higher orders. 



The generalization of (B.5 ) to the case of an arbitrary partial- wave is analogous 



to (B.4), where a subtraction of the following form 



ra^j—lm^ 
rriff—m 

+ 



f(w) dw 
Pcm(w 2 ) (w — m N + m lv ) n w — y/s 

f(w) dw 



1=1 



mjv— 2m, 



Pcm(w 2 ) (y/s -m N + m^Y w - y/s 



Ri 

•s — m N + m/) 1 
(B.6) 



is useful. The function f(w) is regular at the pseudo threshold. The residui 
Ri in (B.6) are to be adjusted such that the r.h.s. of (B.6) does not have any 
contributions from isolated poles at the pseudo threshold y/s = — m^. 



We turn to the u-channel nucleon and t-channel exchange processes in Comp- 
ton scattering. For an s-wave angle average we derive the representations 



i 



dcosO 1 r — 1 ds' 



i 



2 u — m 2 N J 4p 2 m (s') s' — s s — m 2 N ' 



dcosO 1 r 1 dm 2 dC + {m 2 



(B.7) 



2 t-ml J [ 4p 2 m C + (m 2 t ) C + {m 2 t ) - s dm 2 



-1 m% 
oo 



+ 



1 dm 2 dC-(m 2 ) 

4Pc m C-(m 2 ) C_(m?) - s dm 2 



The pole term, singular at y/s — mjv, from the nucleon u-channel exchange 
contributes to Ui ns id e (y/s)- The t-channel one-pion exchange contributes in the 
region with m n < m t < 2m,. Right at m t = 2m n the functions C±(m 2 ) touch 
the contour line C4 of Fig. [3] The generalization to the case of non-s-wave 



type contributions is analogous to (B.4). 



70 



References 



[I] V. Bernard, N. Kaiser, U.-G. Meissner, Chiral dynamics in nucleons and nuclei, 
Int. J. Mod. Phys. E4 (1995) 193-346. 

[2] V. Bernard, Chiral Perturbation Theory and Baryon Properties, Prog. Part. 
Nucl. Phys. 60 (2008) 82-160. 

[3] V. Pascalutsa, M. Vanderhaeghen, S. N. Yang, Electromagnetic excitation of 
the Delta(1232)-resonance, Phys. Rept. 437 (2007) 125-232. 

[4] V. Bernard, N. Kaiser, U.-G. Meissner, Determination of the low-energy 
constants of the next-to- leading order chiral pion nucleon Lagrangian, Nucl. 
Phys. A615 (1997) 483-500. 

[5] N. Fettes, U.-G. Meissner, S. Steininger, Pion nucleon scattering in chiral 
perturbation theory. I: Isospin-symmetric case, Nucl. Phys. A640 (1998) 199- 
234. 

[6] N. Fettes, U.-G. Meissner, Pion nucleon scattering in chiral perturbation theory. 
II: Fourth order calculation, Nucl. Phys. A676 (2000) 311. 

[7] V. Bernard, N. Kaiser, U. G. Meissner, Threshold pion photoproduction in 
chiral perturbation theory, Nucl. Phys. B383 (1992) 442-496. 

[8] V. Bernard, N. Kaiser, U.-G. Meissner, Neutral pion photoproduction off 
nucleons revisited, Z. Phys. C70 (1996) 483-498. 

[9] V. Bernard, N. Kaiser, U.-G. Meissner, Chiral corrections to the Kroll- 
Ruderman theorem, Phys. Lett. B383 (1996) 116-120. 

[10] H. W. Fearing, T. R. Hemmert, R. Lewis, C. Unkmeir, Radiative pion capture 
by a nucleon, Phys. Rev. C62 (2000) 054006. 

[II] S. R. Beane, M. Malheiro, J. A. McGovern, D. R. Phillips, U. van Kolck, 
Compton scattering on the proton, neutron, and deuteron in chiral perturbation 
theory to 0(Q 4 ), Nucl. Phys. A747 (2005) 311-361. 

[12] O. Scholten, A. Y. Korchin, V. Pascalutsa, D. Van Neck, Pion and photon 
induced reactions on the nucleon in a unitary model, Phys. Lett. B384 (1996) 
13-19. 

[13] O. Krehl, C. Hanhart, S. Krewald, J. Speth, What is the structure of the Roper 
resonance?, Phys. Rev. C62 (2000) 025207. 

[14] G. Penner, U. Mosel, Vector meson production and nucleon resonance analysis 
in a coupled-channel approach for energies rriN < sqrt(s) < 2 GeV II: photon- 
induced results, Phys. Rev. C66 (2002) 055212. 

[15] A. M. Gasparyan, J. Haidenbauer, C. Hanhart, J. Speth, Pion nucleon scattering 
in a meson exchange model, Phys. Rev. C68 (2003) 045207. 



71 



[16] V. Shklyar, H. Lenske, U. Mosel, G. Penner, Coupled-channel analysis of the 
omega-meson production in pi N and gamma N reactions for cm energies up to 
2-GeV, Phys. Rev. C71 (2005) 055206. 

[17] A. Matsuyama, T. Sato, T. S. H. Lee, Dynamical coupled-channel model of 
meson production reactions in the nucleon resonance region, Phys. Rept. 439 
(2007) 193-253. 

[18] D. Drechsel, S. S. Kamalov, L. Tiator, Unitary Isobar Model - MAID2007, Eur. 
Phys. J. A34 (2007) 69-97. 

[19] A. Datta, S. Pakvasa, Pion nucleon phase shifts in heavy baryon chiral 
perturbation theory, Phys. Rev. D56 (1997) 4322-4327. 

[20] P. J. Ellis, H.-B. Tang, Pion nucleon scattering in a new approach to chiral 
perturbation theory, Phys. Rev. C57 (1998) 3356-3375. 

[21] N. Fettes, U. G. Meissner, Pion nucleon scattering in an effective chiral field 
theory with explicit spin-3/2 fields, Nucl. Phys. A679 (2001) 629-670. 

[22] M. F. M. Lutz, E. E. Kolomeitsev, Relativistic chiral SU(3) symmetry, large 
N(c) sum rules and meson baryon scattering, Nucl. Phys. A700 (2002) 193-308. 

[23] K. Torikoshi, P. J. Ellis, Low energy pion nucleon scattering in the heavy baryon 
and infrared schemes, Phys. Rev. C67 (2003) 015208. 

[24] V. Pascalutsa, D. R. Phillips, Effective theory of the Delta(1232) in Compton 
scattering off the nucleon, Phys. Rev. C67 (2003) 055202. 

[25] V. Pascalutsa, M. Vanderhaeghen, Chiral effective-field theory in the 
Delta(1232) region. I: Pion electroproduction on the nucleon, Phys. Rev. D73 
(2006) 034003. 

[26] V. Pascalutsa, The Delta(1232) Resonance in Chiral Effective Field Theory, 
Prog. Part. Nucl. Phys. 61 (2008) 27-33. 

[27] B. Long, U. van Kolck, ttN Scattering in the A(1232) Region in an Effective 
Field Theory, Nucl. Phys. A840 (2010) 39-75. 

[28] U.-G. Meissner, J. A. Oiler, Chiral unitary meson baryon dynamics in the 
presence of resonances: Elastic pion nucleon scattering, Nucl. Phys. A673 (2000) 
311-334. 

[29] A. Gomez Nicola, J. R. Pelaez, Unitarized pion nucleon scattering within heavy 
baryon chiral perturbation theory, Phys. Rev. D62 (2000) 017502. 

[30] A. Gomez Nicola, J. Nieves, J. R. Pelaez, E. Ruiz Arriola, Improved unitarized 
heavy baryon chiral perturbation theory for pi N scattering, Phys. Lett. B486 
(2000) 77-85. 

[31] A. G. Nicola, J. Nieves, J. R. Pelaez, E. R. Arriola, Improved unitarized heavy 
baryon chiral perturbation theory for pi N scattering to fourth order, Phys. 
Rev. D69 (2004) 076007. 



72 



[32] S. Weinberg, Nuclear forces from chiral Lagrangians, Phys. Lett. B251 (1990) 
288-292. 

[33] S. Weinberg, Effective chiral Lagrangians for nucleon - pion interactions and 
nuclear forces, Nucl. Phys. B363 (1991) 3-18. 

[34] M. F. M. Lutz, E. E. Kolomeitsev, Baryon resonances in the hadrogenesis 
conjecture, in preparation. 

[35] C. H. M. van Antwerpen, I. R. Afnan, A Gauge invariant unitary theory for 
pion photoproduction, Phys. Rev. C52 (1995) 554-567. 

[36] S. W. MacDowell, Analytic Properties of Partial Amplitudes in Meson-Nucleon 
Scattering, Phys. Rev. 116 (1959) 774-778. 

[37] S. Mandelstam, Determination of the pion - nucleon scattering amplitude from 
dispersion relations and unitarity. General theory, Phys. Rev. 112 (1958) 1344- 
1360. 

[38] J. S. Ball, Application of the Mandelstam Representation to Photoproduction 
of Pions from Nucleons, Phys. Rev. 124 (1961) 2014-2028. 

[39] N. Nakanishi, Proof of partial- wave dispersion relations in perturbation theory, 
Phys. Rev. 126 (3) (1962) 1225-1226. 

[40] G. Sommer, Present state of rigorous analytic properties of scattering 
amplitudes, Fortsch. Phys. 18 (1970) 577-688. 

[41] C. K. Chen, Generalized partial-wave dispersion relation and the n/d method, 
Phys. Rev. D5 (1972) 1464-1471. 

[42] P. W. Johnson, R. L. Warnock, Solution of the unitarity equation with 
overlapping left and right cuts: a tool for study of the s* and similar systems, 
J. Math. Phys. 22 (1981) 385. 

[43] G. F. Chew, S. Mandelstam, Theory of low-energy pion pion interactions, Phys. 
Rev. 119 (1960) 467-477. 

[44] G. Frye, R. L. Warnock, Analysis of partial-wave dispersion relations, Phys. 
Rev. 130 (1) (1963) 478-494. 

[45] A. S. Rinat (Reiner), M. Stingl, Multichannel N/D approach to nuclear 
reactions amongst light nuclei. I. General theory, Annals of Physics 65 (1971) 
141-180. 

[46] P. D. B. Collins and E. J. Squires, Regge Poles in Particle Physics, Julius 
Springer- Verlag, Berlin, 1968. 

[47] M. Stingl, A. S. Rinat, Partial Wave Dispersion Relations for the Nonrelativistic 
Five Nucleon System, Phys. Rev. C10 (1974) 1253. 

[48] L. Castillejo, R. H. Dalitz, F. J. Dyson, Low's scattering equation for the 
charged and neutral scalar theories, Phys. Rev. 101 (1956) 453-458. 



73 



[49] W. R. Frazer, Applications of Conformal Mapping to the Phenomenological 
Representation of Scattering Amplitudes, Phys. Rev. 123 (1961) 2180-2182. 

[50] J. Kennedy, T. D. Spearman, Singularities in partial-wave amplitudes for two 
ingoing and two outgoing particles, Phys. Rev. 126 (4) (1962) 1596-1602. 

[51] R. Koch, A Calculation of Low-Energy pi N Partial Waves Based on Fixed t 
Analyticity, Nucl. Phys. A448 (1986) 707. 

[52] R. A. Arndt, W. J. Briscoe, I. I. Strakovsky, R. L. Workman, Extended Partial- 
Wave Analysis of piN Scattering Data, Phys. Rev. C74 (2006) 045205. 

[53] R. G. E. Timmermans, T. A. Rijken, J. J. de Swart, Determination of the 
charged pion coupling constant from data on the charge exchange reaction anti- 
pp^ anti-n n, Phys. Rev. Lett. 67 (1991) 1074-1077. 

[54] F. Zachariasen, Relativistic model field theory with finite self-masses, Phys. 
Rev. 121 (6) (1961) 1851-1862. 

[55] G. R. Bart, R. L. Warnock, Empirical n/d potential and castillejo-dalitz-dyson 
poles in the pion-nucleon p(ll) state, Phys. Rev. Lett. 22 (1969) 1081-1084. 

[56] B. C. Pearce, I. R. Afnan, THREE-BODY UNITARITY, THE CLOUDY BAG 
MODEL AND THE ROPER RESONANCE, Phys. Rev. C40 (1989) 220-236. 

[57] E. Matsinos, Isospin violation in the pi N system at low energies, Phys. Rev. 
C56 (1997) 3014-3025. 

[58] U.-G. Meissner, U. Raha, A. Rusetsky, Isospin-breaking corrections in the pion 
deuteron scattering length, Phys. Lett. B639 (2006) 478-482. 

[59] H. Krebs, E. Epelbaum, U.-G. Meissner, Nuclear forces with Delta-excitations 
up to next-to-next- to-leading order I: peripheral nucleon-nucleon waves, Eur. 
Phys. J. A32 (2007) 127-137. 

[60] M. C. M. Rentmeester, R. G. E. Timmermans, J. J. de Swart, Determination 
of the chiral coupling constants c(3) and c(4) in new p p and n p partial-wave 
analyses, Phys. Rev. C67 (2003) 044001. 

[61] M. Procura, B. U. Musch, T. Wollenweber, T. R. Hemmert, W. Weise, Nucleon 
mass: From lattice QCD to the chiral limit, Phys. Rev. D73 (2006) 114510. 

[62] P. Buettiker, U.-G. Meissner, Pion nucleon scattering inside the Mandelstam 
triangle, Nucl. Phys. A668 (2000) 97-112. 

[63] R. A. Arndt, W. J. Briscoe, I. I. Strakovsky, R. L. Workman, Analysis of Pion 
Photoproduction Data, Phys. Rev. C66 (2002) 055213. 

[64] K. M. Watson, Some general relations between the photoproduction and 
scattering of pi mesons, Phys. Rev. 95 (1954) 228-236. 

[65] R. Leukel, PhD thesis, Mainz (2001), http://wwwa2.kph.uni-maiz.de/A2/. 

[66] I. Preobrajenski, PhD thesis, Mainz (2001), 

http : / / wwwa2.kph.uni-maiz.de / A2/ . 



74 



[67] A. Shafi, et al., Measurement of inverse pion photoproduction at energies 
spanning the N(1440) resonance, Phys. Rev. C70 (2004) 035204. 

[68] G. Blanpied, et al., iV — > A transition and proton polarizabilities from 
measurements of p( 7 polarized , 7 ), p(j polarized , n ), and p( 7 polarized , 
vr+), Phys. Rev. C64 (2001) 025203. 

[69] J. C. Bergstrom, et al., Measurement of the H-l (gamma, piO) cross-section 
near threshold, Phys. Rev. C53 (1996) 1052-1056. 

[70] A. Schmidt, et al., Test of low-energy theorems for p(gamma(pol.),pi0)p in the 
threshold region, Phys. Rev. Lett. 87 (2001) 232501. 

[71] V. Bernard, N. Kaiser, U.-G. Meissner, Aspects of near threshold neutral pion 
photoproduction off protons, Eur. Phys. J. All (2001) 209-216. 

[72] E. Korkmaz, et al., Measurement of the gamma p ~l pi+ n reaction near 
threshold, Phys. Rev. Lett. 83 (1999) 3609-3612. 

[73] M. A. Kovash, Total cross-sections for pi- p — > gamma n at 10-MeV to 20-MeV, 
PiN Newslett. 12N3 (1997) 51-55. 

[74] V. Bernard, N. Kaiser, U.-G. Meissner, Chiral symmetry and the reaction 7j> — > 
vr Op, Phys. Lett. B378 (1996) 337-341. 

[75] A. Schmidt, PhD thesis, Mainz (2001), http://wwwa2.kph.uni-maiz.de/A2/. 

[76] J. C. Bergstrom, R. Igarashi, J. M. Vogt, Measurement of the H-l (gamma, piO) 
cross-section near threshold. II: Pion angular distributions, Phys. Rev. C55 
(1997) 2016-2023. 

[77] O. Hanstein, D. Drechsel, L. Tiator, A dispersion theoretical approach to the 
threshold amplitudes of pion photoproduction, Phys. Lett. B399 (1997) 13-21. 

[78] S. S. Kamalov, G.-Y. Chen, S.-N. Yang, D. Drechsel, L. Tiator, piO photo- and 
electroproduction at threshold within a dynamical model, Phys. Lett. B522 
(2001) 27-36. 

[79] C. Fernandez-Ramirez, A. M. Bernstein, T. W. Donnelly, The unexpected 
impact of D waves in low-energy neutral pion photoproduction from the proton 
and the extraction of multipoles, Phys. Rev. C80 (2009) 065201. 

[80] D. Hornidge, Proposal for an experiment, MAMI-A2/6-03. 

[81] A. Berstein, Proposal for an experiment, MAMI-A2/10-09. 

[82] V. Olmos de Leon, et al., Low-energy Compton scattering and the 
polarizabilities of the proton, Eur. Phys. J. A10 (2001) 207-215. 

[83] E. L. Hallin, et al., Compton scattering from the proton, Phys. Rev. C48 (1993) 
1497-1507. 

[84] S. Wolf, et al., Compton scattering by the proton using a large-acceptance 
arrangement, Eur. Phys. J. A12 (2001) 231-252. 



75 



[85] D. Babusci, G. Giordano, A. I. L'vov, G. Matone, A. M. Nathan, Low-energy 
Compton scattering of polarized photons on polarized nucleons, Phys. Rev. C58 
(1998) 1013-1041. 

[86] A. I. L'vov, V. A. Petrun'kin, M. Schumacher, Dispersion theory of proton 
Compton scattering in the first and second resonance regions, Phys. Rev. C55 
(1997) 359-377. 

[87] M. Schumacher, Polarizability of the Nucleon and Compton Scattering, Prog. 
Part. Nucl. Phys. 55 (2005) 567-646. 

[88] V. Lensky, V. Pascalutsa, Predictive powers of chiral perturbation theory in 
Compton scattering off protons, Eur. Phys. J. C65 (2010) 195-209. 

[89] V. Baru, et al., Precision calculation of the ir~ deuteron scattering length and 
its impact on threshold pi-N scattering. 

[90] G. F. Chew, M. L. Goldberger, F. E. Low, Y. Nambu, Relativistic dispersion 
relation approach to photomeson production, Phys. Rev. 106 (1957) 1345-1355. 

[91] M. Jacob, G. C. Wick, On the general theory of collisions for particles with 
spin, Ann. Phys. 7 (1959) 404-428. 

[92] W. A. Bardeen, W. K. Tung, Invariant amplitudes for photon processes, Phys. 
Rev. 173 (1968) 1423-1433. 

[93] A. Donnachie, G. Shaw, Photopion Production, Disperions Relations, and 
gamma-rho-pi Coupling, Ann. Phys. 37 (1966) 333-364. 

[94] R. A. Arndt, R. L. Workman, Z. Li, L. D. Roper, Partial wave analysis of pion 
photoproduction, Phys. Rev. C42 (1990) 1853-1863. 

[95] W. Pfeil, H. Rollnik, S. Stankowski, A partial-wave analysis for proton Compton 
scattering in the delta(1232) energy region, Nucl. Phys. B73 (1974) 166-188. 



76 



